@LEVEK 


/;  V ' - 


JOD  yL~% 


DN A 4612F 


CHARACTERISTICS  OF  STRIATION  STRUCTURE 
AND  EFFECTS  ON  SCINTILLATION 


Science  Applications,  Inc. 
1200  Prospect  Street 
La  Jolla,  California  92037 


May  1978 


Final  Report  for  Period  1 October  1974—28  April  1978 


' uu 

fE* 


CONTRACT  No.  DNA  001-75-C-0033 


APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED. 


THIS  WORK  SPONSORED  BY  THE  DEFENSE  NUCLEAR  AGENCY 
UNDER  RDT&E  RMSS  CODE  B322077464  S99QAXHB05404  H2590D. 


D D C 


Prepared  for 
Director 

DEFENSE  NUCLEAR  AGENCY 
Washington,  D.  C.  20305 


/ 1 


Destroy  this  report  when  it  is  no  longer 
needed.  Do  not  return  to  sender. 

PLEASE  NOTIFY  THE  DEFENSE  NUCLEAR  AGENCY, 
ATTN:  TISI , WASHINGTON,  D.C.  20305,  IF 
YOUR  ADDRESS  IS  INCORRECT,  IF  YOU  WISH  TO 
BE  DELETED  FROM  THE  DISTRIBUTION  LIST,  OR 
IF  THE  ADDRESSEE  IS  NO  LONGER  EMPLOYED  BY 
YOUR  ORGANIZATION. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 


I r REPORT  NUMBER 


jiac  READ  INSTRUCTIONS 

WIJt- BEFORE  COMPLETING  FORM 

2 GOVT  ACCESSION  NO.  3 RECIPIENT’S  CATALOG  number"' 


DNA  4612F 


14.  TITLE  (and  Subtitle) 


V CHARACTERISTICS  OF  STRIATION  STRUCTURE  AND 
I EFFECTS  ON  SCINTILLATION, 

_ ___________ 

I David  L. /Sachs- 
Lewis  M./Linson 
j David  C. /Baxter  \ 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Science  Applications,  Inc.  / 

1200  Prospect  Street 

La  Jolla,  California  92037 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


'5.-  TYR*  OF  RERORTi  & PERIOD  COVERED 


Ti  a PERIOD  COVEREC 

Final  .Report  far— Period 

' 1 October  1974  — 28  April  *78 

, -1 

6_  P£_R£Q RM I jdlR crRHP o Re  number  - 

•/  SAI-78-677-LJ » 

8.  CONTRACT  OR  GRANT  NUMBERfs) 

J DNA  OOl-75-C-0033y 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  & WORK  UNIT  NUMBERS 


NWED  Subta*k\ S99QAXHB054-04 


12  REPORT-DATE  __  . , 

Director  //  /'  May  1978  L / 

Defense  Nuclear  Agency  13  numbfr  of  pages 

Washington,  D.C.  20305  104 

4 MONITORING  AGENCY  NAME  & ADDRESS///  different  from  Controlling  Office)  15.  SECURITY  CLASS  (of  this  report) 


/ ■ L 


UNCLASSIFIED 


15a.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  this  Report) 


•WX'  , 


Approved  for  public  release;  distributional  imited.  

j >/  y~£3o  J * 


17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  in  Block  20,  if  different  from  Report) 


10  SUPPLEMENTARY  NOTES 

This  work  sponsored  by  the  Defense  Nuclear  Agency  under  RDT&E  RMSS  Code 
B322077464  S99QAXHB05404  H2590D. 

19.  KEY  WORDS  (Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 


Scinti 1 lation 
Propagation 
Power  Spectra 
Nuclear  Striations 


Striations 
Bi furcation 
Structured  Ionization 


20  ABSTRACT  (Continue  on  reverse  side  if  necessary  and  identify  hv  block  number) 

The  characteristics  of  striations  that  are  expected  to  be  generated  following 
a high-altitude  nuclear  explosion  and  the  scintillation  effects  they  produce 
are  described.  The  shapes  and  density  profiles  of  the  striations  are  para- 
meterized. The  connection  between  the  statistics  of  signal  scintillation  and 
the  parameterized  structural  characteristics  of  the  striations  is  described 
through  the  generation  of  formulas  by  an  approximate  analytic  treatment  of 
the  power-spectral -density  (°SD)  of  single-size  striations.  The  formulas  are  _ 


DD  I JAN^73  1 47 37*  EDITION  OF  I NOV  S5  IS  OBSOLETE 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWfccn  Data  Entered) 


J 


J '4  - 


UNCLASSIFIED 

T 1 5 

k SECURITY  CLASSIFICATION  OF  THIS  PAGEfWien  Data  Entered) 

20.  ^ABSTRACT  (Continued) 

^tested  by  numerical  propagation  calculations  and  found  to  be  adequate  within 
the  uncertainties  of  the  striation  parameters.  The  uncertainties  in  our 
present  understanding  of  the  application  of  the  gradient-drift  mechanism 
towards  a precise  description  of  the  striation  characteristics  and  their 
development  in  time  are  explored.  A methodology  for  calculating  the  time- 
dependent  phase  PSD  resulting  from  the  evolving  striations  is  introduced  and 
a sample  calculation  appropriate  for  late-time  after  burst  is  carried  out 
expl icitly. 


AeCESSIf  N ‘cr 

NTIS  ' 

-Hv.e  Section 

DDC 

v:'f  Section  O 

Wit: 

□ 

JUS  > 

KST'MjMMMY  COOES 

Dirt.  (WAIL  and/or  SPECIft 


UNCLASSIFIED 

L 

SECURITY  CLASSIFICATION  OF  THIS  PAGErWien  Data  Entered) 

\ 

! 

I 


f 

f 

[ 

? 

[• 


(- 

l 

I 

\ 


i 


I 


TABLE  OF  CONTENTS 


Section  Page 

LIST  OF  ILLUSTRATIONS 2 

LIST  OF  TABLES 4 

1 INTRODUCTION 5 

2 STRIATION  STRUCTURE 9 

2.1  INTRODUCTION 9 

2.  2 STRIATION  PROFILES 13 

2. 3 POWER  SPECTRAL  DENSITY 14 

2.  4 POWER  SPECTRAL  DENSITY  OF  STRIATIONS  . 16 

3 SCINTILLATION  DEPENDENCE  ON  THE  PARA- 
METERS OF  SHAPE  19 

3. 1 CIRCULAR  CONTOUR 19 

3. 2 ELLIPTICAL  CONTOUR 33 

3.3  EFFECTS  FOR  A SINGLE  SIZE 43 

4 TIME -DEPENDENT  STRIATION  DEVELOPMENT  . . 47 

4.  1 DISCUSSION  OF  STRIATION  DEVELOPMENT 

IN  HANE 47 

4.  2 DISCUSSION  OF  TIME -DEPENDENT  STRIATION 

DEVELOPMENT 61 

5 BIFURCATION  MODEL  AND  TIME -DEPENDENT 

POWER -SPECTRAL -DENSITIES 73 

5. 1 BIFURCATION  MODEL  APPLIED  TO  HANE  . . 73 

5.  2 TIME-DEPENDENT  PHASE  POWER -SPECTRAL- 

DENSITIES  82 

5.3  SUMMARY 92 

6 LIST  OF  REFERENCES 97 


1 


LIST  OF  ILLUSTRATIONS 


Figure 

1. 

2. 

3. 

4. 

5. 

6. 

7. 

8. 

9. 

10. 

11. 


Schematic  of  electron  density  contours  

Power  Spectral  Density  (PSD)  of  integrated  phase  shift 

for  circular  n contour  striations  (0=1) 

e 

Contours  of  S.  and  Va  for  circular  n contours  0=1) 
and  a=0  ...  

Contours  of  S.  and  Va  for  circular  n contours  (0  -1) 
and  o=0. 9 

Contours  of  S.  and  Va  for  circular  n contours  (3  =1) 
and  «=1. 0 

Power  Spectral  Density  (PSD)  of  integrated  phase  shift 
for  elliptical  ng  contour  striations  (0  ---3)  and  a=0, 1. . . . 

Contours  of  S.  and  Va  for  elliptical  n contours  (0=3) 
and  o=0  

Contours  of  S.  and  Va  for  elliptical  n contours  (0=3) 
and  o=l 

Power  Spectral  Density  (PSD)  of  integrated  phase  shift 
for  elliptical  n contour  striations  (0  =10)  and  a=0, 1. . . . 

Contours  of  S.  and  Va  for  elliptical  n contours  (0=10) 
and  a=0  . . 

Contours  of  S.  and  Va  for  elliptical  n contours  (0=10) 
and  a=l 

Iso -Pedersen  conductivity  contours  for  x = 1 


13.  Iso -Pedersen  conductivity  contours  for  x = 10 

14.  Isodensity  contours  and  striations  for  a one -dimensional 

cloud  with  a periodic  initial  perturbation 


Page 

11 

20 

23 

24 

25 

35 

37 

39 

40 

41 

42 
52 
52 

63 


2 


Isodensity  contours  for  an  initially  circular  plasma 
cloud  with  no  initial  perturbation 


Isodensity  contours  for  an  initially  circular  plasma 
cloud  with  a 5%  initial  perturbation 

Bifurcation  time  as  a function  of  striation  scale  size  . 

18.  The  two  distributions  in  striation  scale  size 

19.  Total  number  of  striations  as  a function  of  time  . . . 

20.  Time  development  of  the  Power-Spectral-Density 

corresponding  to  an  initial  delta  function 

21.  Time  development  of  the  Power -Spectral -Density 

corresponding  to  an  initial  l/a2  distribution 

22.  Comparison  of  the  initial  Power-Spectral-Densities 

at  500  s 

23.  Comparison  of  the  Power -Spectral -Densities  at  1000  s 

24.  Comparison  of  the  Power-Spectral -Densities  at  1500  s 


81 

83 

85 

87 

89 

90 
93 


25.  Comparison  of  the  Power-Spectral-Densities  at  2000  s . 94 


1.  INTRODUCTION 


High  Altitude  Nuclear  Explosions  (HANE)  deposit  large  amounts  of 
ionization  in  the  earth’s  ionosphere  and  above.  At  late  times  the  ionization 
density  is  low  enough  so  that  radio  signals  from  satellites  to  ground  or  air 
receivers  are  not  absorbed.  At  these  times  much  of  the  ionization  is  structured. 
The  ionization  tends  to  form  clumps  (striations)  which  are  nearly  uniform  along 
the  earth's  magnetic  field  for  large  distances  but  break  up  in  the  direction  per- 
pendicular to  the  field  into  smaller  striations.  The  structuring  of  the  ioniza- 
tion causes  random  refraction  and  diffraction  of  the  radio  signal  producing 
scintillation  in  the  signal  phase  and  amplitude  at  the  receiver  as  the  satellite- 
receiver  line-of-sight  moves  across  the  structure. 

The  understanding  of  the  extent  and  character  of  the  scintillation  is 

important  for  the  planning  of  future  satellite  communication  systems.  The 

Defense  Nuclear  Agency  (DNA)  has  been  studying  this  problem  for  many  years. 

There  have  been  many  studies  of  the  effects  of  various  models  of  the  structure 
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on  the  statistics  of  the  signal  scintillation. 

This  report  represents  an  extension  of  these  studies.  A primary  goal 
is  to  attack  the  problem  of  describing  the  scintillation  effects  to  be  expected 
from  HANE.  Most  previous  studies  of  the  development  of  striations  have  not 
concentrated  on  specifying  the  shapes  and  profiles  of  striations  that  might  be 
expected  following  a nuclear  burst.  In  this  report  models  of  the  structure  are 
parameterized  and  an  attempt  is  made  to  relate  these  models  with  the  structure 
to  be  expected  in  HANE.  The  resulting  effects  on  scintillations  are  calculated. 
Furthermore,  the  relation  between  the  model  values  and  the  physical  processes 
that  cause  structuring  is  discussed.  A new  methodology  for  determining  the 
power-spectral-density  (PSD)  of  the  phase  variations  as  a function  of  time  is 
also  discussed  and  an  explicit  example  is  presented. 
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In  principle,  the  procedure  to  follow  is  straightforward.  The  output  of 


the  hydrocodes  determines  the  electron  density,  neutral  density,  and  neutral  winds 
on  a global  scale  (hundreds  to  thousands  of  kilometers).  From  these  parameters  a 
complete  understanding  of  the  physical  processes  which  cause  structuring  would 
provide  the  capability  of  describing  the  development  of  striations  in  the  microscale 
as  a function  of  time  and  space.  Given  the  parameters  that  characterize  the 
irregularities  of  electron  density,  the  power-spectral-density  of  the  variations 
in  phase  as  a function  of  space  and  time  would  be  determined.  That  information 
is  sufficient  to  calculate  various  important  parameters  which  are  a measure  of 
the  effect  on  propagation,  such  as  scintillation  index  or  the  bit-error-rate  on 
communications  links. 

We  do  not  presently  have  sufficient  knowledge  to  carry  out  the  above 
task  with  certainty.  The  weak  link  is  in  providing  an  accurate  description  of  the 
structure.  In  order  to  provide  this  description  three  basic  things  are  needed: 

1)  a specification  of  the  background  conditions;  2)  an  understanding  of  the  basic 
mechanisms  that  lead  to  the  generation  of  structure;  and  3)  detailed  knowledge 
regarding  the  method  in  which  those  basic  mechanisms  operating  under  the 
specified  background  conditions  lead  to  the  generation  of  specific  structures. 

The  various  two-fluid  hydrocodes  and  their  electrostatic  counterparts 
can  be  considered  state-of-the-art  in  regard  to  specifying  the  background  large- 
scale  conditions.  These  various  codes  provide  the  specification  of  the  distribution 
of  the  neutral  particles,  ionization,  and  the  velocity  of  neutral  particles.  V , on 
the  large  scale.  Various  techniques  are  available  for  calculating  the  ion  motion. 
The  ion  velocity,  must  be  consistent  with  the  electric  field  pattern  that  results 
from  the  application  of  the  continuity  of  current  flow  in  the  magnetized  plasma. 
From  these,  the  component  of  the  slip  velocity,  Vg  = Vn  - V.,  perpendicular  to 
the  magnetic  field  can  be  estimated.  The  specification  of  these  time-varying  para- 
meters by  the  output  of  the  hydrocodes  is  necessary  but  not  sufficient  in  order  to 
specify  the  structure. 

In  this  report  we  will  concentrate  on  the  gradient- drift  mechanism  as 
being  the  mechanism  that  leads  to  the  generation  of  striations  during  the  first 
tens  of  minutes  following  a burst.  Since  this  mechanism  was  first  identified 
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nine  years  ago  , it  has  been  the  subject  of  intense  analytical,  experimental, 
and  numerical  work.  However,  our  understanding  of  the  complex  electro- 
dynamics that  governs  development  of  irregularities  is  not  yet  sufficient  to  be 
able  to  specify  the  structure  precisely.  The  hydrocodes  have  been  able  to 
specify  the  regions  of  space  as  a function  of  time  in  which  significant  growth 
of  structure  can  be  expected.  The  microscale  codes,  operated  principally 
by  NRL,  have  provided  some  of  the  characteristics  of  the  development  of  struc- 
ture for  particular  cases.  Observations  of  some  of  the  characteristics  of 
striations  produced  by  the  gradient-drift  mechanism  have  been  obtained  from 
barium  release  experiments  conducted  in  the  ionosphere.  Basically  the 
results  of  current  research  have  been  to  support  the  validity  of  the  general 
concepts  regarding  the  development  of  striations.  However,  many  specific 
questions  for  which  answers  are  required  in  order  to  accurately  specify  the 
structure, have  at  the  moment  only  speculative  theoretical  answers. 

In  order  to  estimate  the  effects  of  striations  on  propagation,  we  have 
accomplished  three  things.  First,  we  have  modeled  characteristics  of  striations 
that  have  an  impact  on  propagation.  These  characteristics  are  the  shape  of  the 
striation,  the  profile  of  the  electron  density  within  the  striation,  and,  in  some 
cases,  a distribution  of  scale  sizes  of  striations.  Second,  we  have  identified 
specific  gaps  in  our  current  understanding  of  the  physical  processes  that  specify 
these  striation  characteristics.  We  indicate  that  there  are  procedures  that 
could  be  followed  that  could  verify  or  deny  some  of  the  theoretical  speculations 
regarding  the  specifications  of  these  characteristics.  Thirdly,  we  have  developed 
a methodology  for  determining  the  time -dependent  power-spectral-density  of 
striations  and  have  carried  out  a specific  model  calculation  exemplifying  the 
results  of  applying  that  methodology.  This  methodology  can  be  viewed  as  an 
alternative  to  the  methodology  being  developed  by  J.  Workman  of  BRA  for  the 
SCENARIO  code. 

Section  2 contains  a description  of  the  model  structures  used  in  order 
to  calculate  the  effects  of  striations  on  propagation.  The  parameterization  of 
the  shapes  and  profiles  are  described  along  with  the  resultant  power-spectral- 
densities.  The  significance  of  the  power-spectral-density  (PSD)  for  estimating 


the  signal  scintillation  statistics  is  emphasized.  Section  3 contains  the  results  of 
applying  the  various  models  in  order  to  obtain  the  dependence  of  the  signal  scintil- 
lation statistics  on  the  assumed  structures. 


Section  4 discusses  the  physical  processes  that  lead  to  a specification 
of  the  structures.  This  discussion  includes  an  application  of  our  present  under- 
standing derived  from  both  linear  and  nonlinear  analyses  to  a description  of  the 
time -development  of  striations  in  HANE.  The  gaps  in  our  knowledge  for  making 
precise  statements  are  emphasized.  A methodology  based  on  the  bifurcation  of 
striations  for  determining  the  time -dependent  power -spectral -density  is  described 
in  Section  5.  A specific  example  appropriate  for  HANE  conditions  is  calculated 
in  detail.  It  is  emphasized  that  the  inputs  to  this  methodology  are  the  physical 
processes  that  lead  to  the  development  of  striations.  As  our  understanding  regarding 
the  development  of  striations  is  improved,  those  aspects  that  have  a direct  impact  on 
scintillation  effects  can  be  incorporated  into  the  same  methodology. 


ku ' A 


2.  STRIA TIQN  STRUCTIRE 


2. 1 INTRODUCTION 

A radio  signal  passing  through  a medium  will  have  spatial  variation 
in  its  amplitude  (scintillation)  over  an  area  perpendicular  to  the  line-of-sight 
if  the  medium  is  structured,  i.  e. , not  homogeneous.  A late-time  high-altitude 
nuclear  burst’s  ionization  is  structured  in  the  following  ways  on  a descending 
order  of  scale  size: 

1.  The  ionization  is  localized  within  a region  hundreds  of  kilometers 
in  radius  transverse  to  the  magnetic  field  and  thousands  of  kilo- 
meters in  extent  along  the  magnetic  field; 

2.  The  ionization  is  in  the  form  of  distinct  cylindrical  shapes 
(striations)  imbedded  in  a region  of  much  smaller  (background) 
ionization; 

3.  The  boundaries  of  these  striations  (electron  density  contours)  are 
not  smooth; 

4.  Sound-wave -like  disturbances  can  exist  within  the  striations. 

The  effects  caused  by  the  large-scale  structure  of  the  first  type  alcne 
are  negligible  for  the  late -time  electron  densities  and  radio  frequencies  of 
interest  in  this  study. 

The  dependence  of  the  system  degradation  on  the  signal  scintillation 
caused  by  structure  of  the  second  type  is  the  subject  of  concern.  As  in  our  pre- 
vious studies*  we  model  this  structure  as  a given  set  of  ionization  profiles 
randomly  distributed  within  a given  area.  The  amount  of  structure  of  tjpes 
3 and  4 are  unknown  at  this  time.  The  statistics  of  the  scintillating  signal  at 
some  distance  beyond  the  striations  are  obtained  using  the  thin-screen  approxi- 
mation. The  results  are  presented  as  a function  of  A,  the  signal  wavelength, 
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d,  distance  away  from  the  striations,  nQ,  peak  electron  density,  and  number,  size, 
and  shape  of  the  ionization  profile. 

For  any  given  set  of  striation  sizes  and  shapes  the  signal  and  its  statis- 
tics are  functions  of  only  two  dimensionless  parameters: 


where  a is  the  characteristic  scale  size  of  the  striations;  and 

<Pq  =J(tf>(x)  - 0)2 

the  root-mean-square  of  the  phase  variations  in  the  signal  at  the  exit  plane  of  the 
striations.  In  the  thin-screen  approximation,  this  phase  is  obtained  by  a straight- 
line  integration  of  the  electron  density  along  the  line-of-sight  of  the  radio  signal. 
For  the  electron  density  defined  in  an  area  L by  L in  the  x-y  plane; 


fL 

<p(x)«r  X J dyn  (x,y)  . 


where  y is  the  direction  of  propagation  of  a plane  wave  of  wavelength  \ and 

r = 2.  818  x 10' 13  cm 
e 

is  the  classical  electron  radius.  We  define  the  average  of  any 
function  f(x)  as 


dx  f (x) . 


We  take  the  x-direction  as  the  direction  of  motion  of  the  line-of-sight  (y-direction) 

of  the  signal  past  the  ionization  due  to  the  motion  of  the  transmitter,  receiver  or 

2 

striations.  The  area  A = L , will  be  expanding  in  time  as  the  ionization  spreads 
across  magnetic  field  lines. 

Consider  the  following  typical  example  of  the  time-dependent  ionization 
structure  in  the  plane  of  area  A.  Before  the  onset  of  striations,  assume  that 
the  ionization  has  a radial  dependence  only  about  some  central  position.  Fig.  1(A) 
illustrates  the  contours  of  constant  electron  density,  ng,  for  this  case.  In  the 
general  case,  the  magnetic  field  will  not  be  perpendicular  to  the  plane  and  the 
contours  will  be  ellipses  elongated  in  the  direction  of  the  magnetic  field’s  pro- 
jection on  the  plane.  This  complication  does  not  qualitatively  change  the  result. 
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The  contours  illustrated  are  a representation  of  the  output  of  an  MHD 
computer  code  minutes  after  a high  altitude  burst.  These  contours  can  be 
approximately  fit  by  a function  of  the  form 


2,  2 

ne(x,  y)  = nQ  e"r  ao 


+ n, 


where  r is  the  distance  from  the  center  and  aQ  » 150  km  for  a typical  case. 

This  large  value  of  a characteristic  size  aQ,  does  not  lead  to  any  ampli- 
tude scintillation  in  a signal  traversing  the  ionization  for  the  range  of  parameters 
of  concern  to  this  study: 

X s 10~3  km 

n £ 103  cm  3 

n.  £ 10°  cm*J 

b 

We  therefore  ignore  this  structure  of  the  first  type. 

When  the  ionization  breaks  up  from  instabilities  generated  by  a radial 
neutral  wind  into  fingers  that  are  convected  radially,  the  total  electron  content 
will  be  maintained  except  for  the  ongoing  processes  of  chemistry  and  diffusion 
along  the  magnetic  field.  The  fact  that  the  ionization  is  structured  into  smaller 
parcels  will  bring  about  amplitude  variations.  This  is  the  structure  of  the  second 
type  to  be  investigated  in  this  report. 

Barring  unforeseen  surprises  we  expect  the  late-time  nuclear  structure 
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to  resemble  a radial  version  of  the  output  of  the  electrostatic  codes  of  NRL  . By 
radial  we  mean  that  the  2-D  cartesian  fingers  of  NRL  which  extend  in  the  constant 
neutral  wind  direction  would  appear  in  the  nuclear  case  to  extend  radially  out- 
ward from  a cylinder  axis  along  the  magnetic  field  through  the  center  of  the  ioni- 
zation [Fig.  1(B)].  in  the  course  of  time  there  will  be  bifurcations  which  increase 
the  number  of  fingers  and  decrease  their  size.  The  orientation  of  the  major  axes 
of  the  fingers  will  be  primarily  along  radials  and  at  some  point  further  bifurcation 
will  cease  and  there  will  be  a limiting  size  distribution  with  the  smallest  being 
depleted  by  diffusion.  A model  of  this  time -dependent  process  will  be  described 
in  more  detail  in  Section  5. 


i 
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STRIA TIGN  PROFILES 


The  essential  properties  of  the  structured  ionization  of  a single  striation 
important  for  producing  amplitude  scintillation  appear  to  be: 

1.  Shapes  which  vary  from  circular  to  highly  elongated  ellipses; 

2.  Profiles  with  large  density  gradients  in  the  narrow  dimension; 

3.  Smaller  gradients  in  the  elongated  direction  down  to  some  value 
of  ng,  then  abrupt  gradients. 

We  choose  to  model  these  features  with  a two-parameter  set  of  analy- 
tic profiles: 


/ 2 2. 

"e  ~ "o 


r s a 


= 0 r > a (2.2.1) 

2 2 2 

where  r = 0x  + y / /3 . Contours  of  electron  density  are  ellipses  with  major  to 
minor  axes  ratio  p . 

This  is  a new  concept  of  profile  which  occurred  to  the  authors  after 

g 

seeing  the  NRL  contours  of  n from  their  numerical  codes  . In  contrast  to  earlier 
7 e 

work  , this  numerical  study  treated  a case  where  the  initial  Pedersen  conductivity 
was  much  higher  than  the  background.  The  profiles  of  ionization  varied  smoothly 
down  from  the  peak  value  to  some  smaller  value  where  the  fall-off  became  abrupt. 
Physical  motivations  for  this  type  of  behavior  will  be  discussed  in  Section  4. 

We  use  the  parameter  a to  characterize  the  value  of  density  at  which 
the  profile  becomes  steep.  The  parameter  p describes  the  ellipticity  of  the 
density  contours.  The  parameter  a is  the  characteristic  size  of  the  profile.  If 
the  striations  are  distributed  over  a large  range  of  sizes  a,  the  form  of  the 
distribution  of  sizes  will  be  another  parameter  of  the  problem.  We  have  been 
unable  to  determine  from  the  phenomenology  of  the  bifurcation  of  striations  any 
particular  form  for  the  distribution.  We  have  shown  in  previous  work*  that  the 
scintillations  are  highly  dependent  on  the  profile  shape.  A possible  expected 
distribution  of  sizes  will  be  discussed  in  Section  5.  The  comparison  of  the 
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distribution  of  sizes  and  the  single  size  will  be  made  by  means  of  the  power- 
spectral-density  of  the  electron  density. 

2. 3 POWERrSPECTRAL-DENSITY 


We  have  chosen  the  basic  representation  of  the  structured  electron  density 
as  a set  of  striations  of  given  number,  size  and  profile.  An  alternate  representation 
of  the  electron  density  is  through  its  power- spectral -density  (PSD). 

2 

From  the  Fourier  transforms  of  the  electron  density  in  the  area  A = L , 


f©  = /L  */Lf 


e-ik‘xne(x,y)  , 


we  have 


where 


n (x,y)  = £ f(k)elk'X 
e 1,  m 


; = (■£!_.  *5L) 


» <1,  m <» 


The  average  electron  density  is 

ne  = y J dx  ng(x,  y)  = f(0)  . 


If  we  set  f(0)  a 0,  then  f(k)  is  the  Fourier  transform  of  An  * n - n . Define 

” C v 

the  PSD  of  the  electron  density  by 


*(S)  =-^  y**c  e"lk*  x Ane(x’ 1 A ne(x’  + x)  (2.3.1) 

= \ ft f<&'  e'lk  X Ang(x’ )Ang(x  ’ +x) 

A 

= f(E)i(-E)  = I f (C>  j 2 

where  f(-k)  = f*(k)  follows  from  the  fact  that  Ang  is  real. 
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Thus  the  PSD  of  ng,  4(k),  defined  by  Eq.  (2. 3. 1),  is  given  by 

2 

*(k)=|f(k)l  . 

The  propagation  analysis  requires 

rh 

«(x)  E r Xj  An(x,  y)dy 


which  has  the  Fourier  transform 

F(k)  - i J 


iloc 

<p(x)e  dx  = rg  XL  f ( k,  0)  . 


8 


We  will  use  the  notation  of  Salpeter  for  the  PSD  of  phase 

00 

cp2(q)  ■ f dre‘iqr  Q(x)Q(x  + r) 


defined  for  an  infinite  interval  and  evaluate  it  for  our  defined  phase  on  the 
interval  0 s x s L to  obtain 

o2(q)  | F(q)!2  (reXL)2  ]f(q.  0)  |2  (2.3.2) 

8 9 

In  Salpeter’ s description  of  the  original  analysis  of  Mercier  , the 

statistics  of  the  intensity  scintillations  are  described  completely  in  terms  of 
2 

to  (q).  This  is  because  Mercier  required  the  condition 

ei  Ea.  0(x.)  = e-l[Ea.  <fl(x.)]2 
1 ] 

for  a.  any  set  of  constants.  This  condition  is  satisfied  if  <p(x.)  is  a Gaussian 
1 3 

distributed  variable.  Much  of  the  recent  literature  has  been  devoted  to  the 

2 

question  of  the  requirements  for  <p  (q)  to  be  sufficient  for  the  description  of 

signal  scintillations.  A sufficient  requirement  is  that  the  phase  be  the  sum  of 

many  random  contributions.  Since  this  is  the  case  in  our  analysis  we  will 

2 

consider  that  a specification  of  <p  (q)  is  sufficient. 


r 


I 


m 
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2.  4 POWER-SPECTRAL -DENSITY  OF  STRIA TIONS 


follows. 


The  connection  between  <p  (q)  and  a random  set  of  striations  is  as 
The  profile  in  expression  (2. 2. 1)  has  the  Fourier  transform 

— » -4 

-ik  • X; 


f(kx„  ky,;  ;,)  = e-  • if,(kx,,ky,) 


where 


2 

2n  if  a f J i (s)  2 2 a 2 

«,V  V-ff-i-T-  ) + 


J0(S) 


where 


/‘V  2 

s = * J-j-  +V* 

and  J.  0 are  Bessel  functions. 

■4 

The  y’-direction  is  along  the  major  axis  of  the  ellipse  and  is  the 
position  of  the  striation  in  the  xy  plane.  The  phase  integral  is  along  the  y-axis 
which  is  taken  to  be  at  an  angle  0 ^ to  the  y'-axis  so  that 

x = x’  cos  0^  + y’  sin  0.  . 


The  Fourier  transform  of  the  phase  produced  by  the  striation  is 
FjOO  - reXLe-‘tal^S  f e-ikxns(x-x..  y-y,) 

= rciLe-iltxi  e "i(k  C0S  V + k Sto  y'-y.) 


= rp>Le’ikXi  f.(k  cos  0.,  k sin  0.)  . 
c 111 

The  total  phase  produced  by  a set  of  striations  has  the  Fourier 


transform 


F(k)  = E F.(k) 
i 


16 


and  the  PSD 


<p  (q)  --2^(re*^S  f^  (q  cos  0.,,  q sin  0.) 

+ s e“iq(xi“xP  f^q  cos  0.,  q sin  0.)  f.(q  cos0.,  q sin  0.) 

1 x j ] ] j 


(2.  4. 1) 


from  Eq.  (2. 3.  2). 


. 1 


Some  comments  are  in  order  here.  In  a previous  analysis  we  showed 
the  contribution  of  the  second  term  in  brackets  to  be  an  A.  C.  jitter  about  the 


contribution  of  the  first  term  due  to  the  positions  of  a given  set  of  striations.  If 
2 

we  average  <p  (q)  over  an  ensemble  of  random  placements  of  the  striations  with 


the  fractional  number  of  times  striation  i is  at  xi  being  P^x.Jdx.  and  consider 


the  positioning  of  striations  to  be  uncorrelated,  the  average  of  the  second  term  in 
brackets  is 


s fi(q  cos  ft,  q sin  0.)f{(q  cos  0. 
tfj  1 1 ] ] 


, q sin  0 ^ IS 


where 


= f P^x.) dx.  e'iqxi  . 


If  the  distribution  over  the  length  L is  uniform,  P = - and  S = ^^fak/2) 


which  is  zero  at  q = 2*  — t = 1,  2,  . . . , the  only  values  of  q used. 

i-i 

This  means  the  PSD  is  given  by  the  first  term  alone  for  determining 
effects  due  to  the  striations'  profiles. 


There  have  been  experimental  determinations  of  the  PSD  which  have 
included  large  contributions  from  structure  of  the  first  type,  the  positioning  of 
the  ionization  in  space.  These  contributions  can  be  illustrated  by  distributing 


the  striations  over  the  line  L by  a Gaussian  probability  P(x.)  where 


-(x.-L/2)2/a2 

P(x.)  = — a « L/2  . 

1 J¥T 
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j 

J 


1 


2 2/4 

Then ! S | = e~q  a ' . if  there  are  M striations  with  the  same  form  of  profile, 

the  bracket  in  Eq.  (2. 4. 1)  becomes 

2 2 

Mf2  [1  + (M-l)e-q  a ^2] 

- 2 2 /2 

and  the  second  term  dominates  at  low  values  of  q until  Me"q  a 6 becomes  less 
than  unity.  This  is  a real  effect  which  must  be  included  if  one  wants  to  include 
the  effects  of  structure  of  the  first  type  on  scintillation  or  if  one  wants  to  extract 
properties  of  striations  from  experimentally  obtained  PSDs. 


Another  point  to  be  made  is  that  the  portion  of  the  PSD  due  to  structure 
of  the  first  type  does  not  satisfy  the  criterion  for  the  PSD  to  be  sufficient  to 
describe  scintillation  effects.  That  is,  the  phase  variations  resulting  from  this 
type  of  structure  are  not  the  resultant  of  many  random  inputs. 


For  our  purposes,  the  connection  between  striation  profiles  and  the 
PSD  will  be  through  the  first  term  alone.  Then 

o2(q)  =^~  (reX)2  E f2 (q  cos  0.,  q sin  0.)  . 
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SCINTILLATION  DEPENDENCE  ON  THE  PARAMETERS  OF  SHAPE 
CIRCULAR  CONTOUR  (ft  = 1) 

We  first  discuss  the  case  of  fi  = 1 (circular  ng  contours).  The  profile  is 


/ 2 2 . 


r £ a 


where  a = 0 produces  the  rod  profile  studied  previously  and  a = 1 produces  a 
parabolic  profile.  The  PSD  for  this  profile  is 


2,  . 2ff  / 2,2  rJl^  M 2,  „ 2 J2^1 

^ (q)  = T"  (reXV  > L“q ~ (1  " “ } + 2a  j 


(3.  1.  1) 


Fig.  2 shows  the  PSD  of  three  profiles  corresponding  to  a = 0,  0.  9 and  1.  After 
the  first  zero  of  the  Bessel  function,  the  locus  of  peaks  is  plotted. 

-3 

The  large  q behavior  of  a = 0 and  0.  9 is  q where  the  large  q behavior 
-5 

of  a = 1 is  q . The  PSD's  are  normalized  so  that  the  profiles  have  the  same  nQ 
and  the  same  total  electron  content  which  requires 

tra^(l  - <x  /2) 

to  have  the  same  value  for  any  o or 

a = a (0)  (l-aV2)_*. 

They  therefore  have  the  same  value  as  q - 0.  The  a = 0.  9 curve  looks  the  same 


as  the  a = 1.  0 curve  until  a region  around  the  value  of  qa  (0)  = 2d 


2</ 1 - a_/2 
: (1  - «) 


after  which  the  q behavior  takes  over.  The  curves  are  all  plotted  vs. 
v = for  the  case  a(0)  = 0.  282.  so  the  break  point  is  at  v = ^57^2)  = 

From  the  curve  we  see  the  transition  to  begin  at  the  second  peak  of  the  Bessel 

_4 

function  where  the  line  has  a v dependence.  The  transition  ends  around  v = 7. 4 
-3 

where  the  u dependence  takes  over. 
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Our  analysis  of  the  scintillation  effects  is  through  the  thin  screen  approxi- 
mation and  a Fast  Fourier  Transform  propagation  routine  as  described  in  our  last 
analysis*.  In  that  analysis  we  obtained  the  required  p(x)  by  summing  contributions 
from  randomly  placed  striations.  In  the  present  analysis  we  have  compared  the 
results  of  that  procedure  with  the  alternate  method  of  using  the  PSD  directly.  The 
alternate  method  stipulates  the  Fourier  transform  of  p(x)  to  be 

• «■)  e19(<l), 

where  the  amplitude  is  the  square  root  of  the  PSD  and  the  phase,  8(q),  is  chosen 
randomly  between  0 and  2ff  for  each  value  of  q.  The  two  methods  agree  if  there 
are  enough  randomly  placed  striations  so  that  the  generated  p(x)  is  due  to  about 
ten  or  more  striations  contributing  for  each  value  of  x.  If  there  are  fewer  stri- 
ations of  larger  electron  density  so  that  p is  unchanged,  the  comparison  shows 

rms 

weaker  scintillation  effects  using  the  random  striations  than  using  the  PSD.  This 
follows  from  the  fact  that  the  PSD  method  assumes  the  phase  to  be  distributed  as 
randomly  as  possible  with  the  same  p g and  correlation  lengths.  The  PSD 
method  is  consistent  with  our  expectation  that  signals  will  pass  through  many 
random  striations. 


Since  the  PSD  is  an  even  function  of  q,  its  Fourier  transform,  the  phase 
autocorrelation,  is 

. . . . r 2 2 

p(x)p(x+r)  = 2 J dqp  (q)  cos  (qr)  sp  p(r) 


where  the  variance 


p2  = 2 

o 


“ « 16  (r  Xn  a ) 

f dq„2(,)=  6 0 


I",,  2V2  16a2(l-a2)  320,’  1 

[(l-o  ) + 


lair 


The  dependence  of  the  phase  variance  on  the  parameter  a is  slight.  We 
maintain  the  electron  content  constant  as  we  vary  a so 

2,  vf  , . _2,2  2 2 


<P  (q=°)  =-217  (reXn0a  ) (l-«  / 2) 


is  constant  and 
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2 _i  2 

Using  a(o)  = a(0)(l-a  /2)  , the  dependence  of  « on  a alone  varies  from  a 

32 

factor  of  unity  for  o = 0 to  a factor  r =0.  862  for  <*=  1.  This  dependence 

105  (i) 2 

*r 

of  <p  =V(fl  o on  a will  be  ignored  and  we  will  present  our  results  of  the  pro- 
pagation study  in  the  usual  manner1  as  a function  of  the  two  parameters,  <po  « o g 

and  e = — M — • 

4 ira^(0) 

One  measure  of  the  severity  of  signal  intensity  scintillation  is  the 
scintillation  index  defined  by 

S?  ■?-  1 

4 

where  I is  the  signal  intensity  on  a plane  a distance  d from  the  striation  normalized 
to  intensity  in  the  absence  of  striations.  The  average  intensity  1 = 1 by  conserva- 
tion of  energy.  This  means  no  energy  is  backscattered  by  the  striations  using 
our  propagation  methods  which  are  strictly  valid  for  small-angle  forward  scatter. 
We  have  found  in  a previous  study1  that  for  S4  s 0.3,  the  scintillation  is  too  mild 
to  disturb  communication  systems  and  for  * 0.9  severe  scintillation  occurs 
which  can  be  approximated  by  Rayleigh  statistics.  These  values  of  are  there- 
fore useful  for  the  discussion  of  system  effects.  Another  useful  parameter  is 
the  intensity  correlation  length,  4.  defined  where 

l(x)l(x+  f)  - 1 = e-1S42  . (3.1.2) 


We  plot  the  contours  of  S4  = 0. 3 and  0.  9 on  Figs.  3,  4 and  5 for  the  cases  a=0, 
0.  9 and  1 respectively  as  a function  of  e and  <pQ.  We  also  plot  contours  of 
t/a(«)  as  a function  of  fand  oQ  in  the  figures  as  dashed  lines.  We  are  strictly 
interested  in  l only  in  the  regime  of  strong  scintillation  S4  * 0. 9.  From  the 
figures  it  is  obvious  that  in  this  case  t/a(a)  is  well  approximated  as  a function 
of  <0Q  alone.  In  fact  4=  fa (a)/<P0  fits  very  well. 
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(radians) 


Figure  5. 


Contours  of  and  t/a  for  circular  n 
and  a=  1.  0. 


contours  (p  = 1) 
e 


\ 


i 


-jar 
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Mercier's  analysis  of  scintillation  effects  in  the  limit  as  e-«®  obtained 
an  expression  for  I(x)I(x+t) which  depended  on  the  PSD  of  phase  alone.  Again, 
this  is  because  he  assumed  the  phase  to  be  Gaussian  distributed  as  we  do.  His 
expression  is 

2 2 

I(x)I(x+r)-  1 = e'2<po  [e2<po  P^  . l] 

which  results  in 

S2  = 1 - e'2tp°2 
4 

or  S.  =0. 3 for  =0.  217  and  S.  =0.  9 for  o =0.  91.  These  values  are  in  agreement 
4 o 4 o 

with  the  top  of  the  contours  in  Fig.  3 for  a = 0.  For  the  other  figures,  the 
S^  = 0. 3 contour  is  in  agreement  but  the  = 0.  9 contour  is  at  a larger  value  of 
0Q  but  appears  to  be  approaching  Mercier’s  limiting  value. 

The  curvature  of  the  contour  of  as  they  approach  Mercier’s  c = ® 

limit  contrasts  strongly  with  the  contours  of  i/a  which  are  straight  lines  down  to 

the  S.  = 0.  9 contour  even  when  * « 1.  We  therefore  believe  Mercier's 
4 

analysis  to  be  valid  for  S4  v 0.  9 for  any  e.  Since  this  regim^  always  has 

0 >1,  we  take  the  limit  of  Mercier's  expression  where  e~2<^o  « 1 and  obtain 
° 2 2 

I(x)I(x+i)  - l=e-2K>  " *0  From  our  formula  (3.  1.2)  for  twe  have 

2 2 1 
<e0  -0O  p(*)  = r - 

2 

Since  we  restrict  ourselves  to  the  case  of  <pn  » 1 we  need  to  evaluate  p(i)  for 

2 2 U 

small  values  of  i/a  where  p * 1 - Ai  /a  where  A is  some  factor  depending  on  the 
functional  form  of  the  auto  correlation  p(i).  Thus 

f ^ a _ 1 _a 

" v^A"  % 

and  the  dependence  of  i/a  on  «oo  is  as  seen  in  the  figures. 

We  shall  now  discuss  the  S^  = 0.  3 contour.  All  figures  show  the  con- 
tour to  asymptote  to  a straight  line  0Q  = constant  for  c » 1 and  to  asymptote  to 
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a straight  line  e = B<pQ~c  for  e « 1 which  always  is  in  the  regime  oQ  »1.  The 
case  a = 0.  9 shows  a break  at  very  small  e from  the  straight  line. 

9 

This  behavior  is  in  agreement  with  the  analysis  of  Mercier  who 


showed  that  for  o « 1, 
o 


In  the  lim  S„ 


00 

= 8 f dq  cp2(q)  sin2  (fq2a2)  . 


/■ 

2 2 

dq  <p  (q)  = 2toQ  or  S4  - J2<P0  • 


It  has  since  been  shown  by  Jokipii10  and  others  that  Mercier' s 

2 

formula  is  also  valid  when  oq  > 1 so  long  as  the  resultant  value  of  « 1.  The 
small  e limit  produces 


~2  2 o4 

= o(  a 


f dq  o2  (q)  q4 


which  diverges  for  the  power  spectra  we  are  considering.  If  the  power  spectra 

-5 

fell  off  faster  than  q for  some  large  value  of  q as  they  must  from  diffusion, 

there  would  be  some  small  ? below  which  the  curve  would  asymptote  to  the 

straight  line  c = B/<pq.  We  assume  this  value  of  e to  correspond  to  a diffusion 

cut  off  size  an  ~ 30  m.  We  shall  assume  our  curves  of  the  power  spectrum  are 
u _1 

accurate  out  to  q = (30  m)  . It  is  interesting  that  the  expansion  of 

2 2 2 2 
sin  (cq  a ) as  a means  of  obtaining  the  small  e limit  which  leads  to  a rapid  e 

2 

fall  off  for  (as  was  found  for  an  exponential  power  spectrum  in  a previous 
study  *)  may  be  incorrect  for  power  law  spectra  even  if  a reasonable  diffusion 
cut-off  exists. 

2 

We  shall  obtain  the  small  c behavior  of  S.  by  assuming  that  the 
2 2 2 4 

sin  (t  q a ) term  wipes  out  the  contribution  to  the  integral  of  the  low  qa 

2 2 2 

portion  of  <p  (q).  (valid  when  cq  a « 1)  and  that  we  can  go  directly  to  the 
asymptotic  limit  of  the  Bessel  functions,  where  for  a = 0 we  have 


„ r . . 2 . 2 2.  4 . 2. 

g2  _ 8 J dq  sin  (eq  a ) y (r^a  ) 


2. 2 sin2  (qa-  ff/4) 


00 

2 f dq 

■ <fl  / — 

0 / \3 

O (qa) 


sin2  (cq2a2)  sin2  (qa-ff/4)  . 


Let  t = cqV.  Then 


where  for  small  e , 


S4 ' 3 f »o2  / dt  (-^r1-)2  sl"2  - »/« ) 
sin2  (yt/7-  ff/4  y 


through  the  important  part  of  the  integrand  so 


00  9 

„2  3 2 / ,,  (slTitY  3ft  2 

s4  = 5 e”0  J dt\— ; = t f®o  • 


2 2 

Therefore  S.  falls  off  as  c rather  than  e . If  there  were  an  abrupt  diffusion  cut- 

2 ’ 2 2 
off  of  <0  (q)  at  q = q then  the  result  is  unchanged  if  eq  a > 1 as  mentioned  above. 

2 2 u 2 U 
For  e<lD  a < 1 we  recover  the  e behavior. 

The  = 0.  3 contour  in  Fig.  3 agrees  very  well  with  a combination  of 
the  two  limits  like 


s2  2B° 
s4 " z s"T 

{l  + Wl] 


(3.  1.3) 


for  example.  The  = 0. 9 contour  in  Fig.  3 has  a high  e limit  which  corresponds 
to  Mercier’s  = 1 - e"2(po  formula.  If  we  arbitrarily  extend  the  formula  of 
Mercier  to  all  c we  obtain 


S2  = 1 -e"2<0O2/(1+ 
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2 

The  = 0. 9 contour  from  this  formula  corresponds  to  = 1.  66  in  Eq. 

(3. 1.3 ).  So  if  we  would  plot  <p=  0.  9 J 1 + for  the  S = 0.  9 contour  and 

€ 4 

<pQ  = 0.  2 J 1 + for  the  = 0. 3 contour  they  would  lie  along  the  plotted 
contours  in  Fig.  3 well  within  the  accuracy  of  the  plots.  Therefore  we  may  say 
for  the  case  of  a = 0,  /J  = 1,  that: 

8~ 

+ 3?7 


2.  scintillation  is  Rayleigh  for  tpQ  > 0.  9 
with  a correlation  length 

l = 0. 6 a/<0Q 

3.  If  desired,  the  amplitude  distribution  for  intermediate  values 
of  could  be  approximated  by  a Ricean  using 

®eff  = V^1+^ 

if  no  attempt  is  made  to  correlate  amplitude  and  phase. 

Our  propagation  output  includes  a calculation  of  the  average  bit  error 
probability,  P , for  noncoherent  frequency  shift  key  (FSK)  modulation  according 
to  the  formula 

p"e  = f P(s)Pe(s)ds 


J7T* 


1.  scintillation  is  negligible  for  < 0. 2 /I 


where  P(s)  is  the  probability  density  that  the  signal -to-noise  ratio  is  s and 
_ «/o 

Pe(s)  = |e~  ' is  the  bit-error  probability  for  the  value  s.  This  is  evaluated 
by  averaging  Pe(Iso)  for  our  output  considering  the  signal-to-noise  ratio  to  be 
sQ  in  the  absence  of  scintillations.  For  (15  dB),  a Rayleigh  distribution 

produces  P"  = . 03  whereas  no  scintillation  means  p"  = 7. 10  ® (or  no 
error).  Our  results  show  Pe  to  vary  with  S^,  e.  g. . a contour  of  Pg  = 10 
follows  the  = 0.  8 contour  very  closely. 
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We  now  skip  to  the  case  a = 1.  Fig.  5 indicates  that  the  small  e behavior 
of  the  S4  = 0. 3 contour  approaches 

S2  = C o 2 e2 
4 0 

whereas  the  = 0. 9 contour  shows  a weaker  e-dependence  and  indicates  a tendency 

to  close  on  the  S . = 0. 3 contour  as  e decreases.  This  behavior  is  similar  to  the 
4 1 

case  of  the  exponential  power  spectrum  studied  earlier  wherein  there  was  focusing 
leading  to  values  of  greater  than  unity. 

For  our  purposes  the  important  aspect  of  the  a = 1 case  is  the  fact  that 
the  = 0. 3 contour  requires  a larger  than  the  a = 0 case,  that  is  <p^  is  larger 
when  c < 1.  For  a=  1 our  approximate  formulas  for  the  small  e case  diverge.  We 
need  to  evaluate 


g2  = 64ir  (r^r^a2)  J dq  -2  -g—  sin2  (cq2a2) 


315* 

Ten 


QO 

2 \ ( dt  . 2 _ 2 

°c  ' J T*sm  (t),2 


O t' 


(jVz)  ■ 


The  integral  can  be  split  into  two  parts  for  small  e 


f'  dt  ti  yr 

J dtt  —i — “ w 


$ sin2  (t)  cos2  (jtJT-  v/i) 


J <3- 2 me) 

e t 
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producing 


_2  2 2 (315)  /.  . 1 v 3 v 

s4  = ”o  f -ir-(ln(T)  + s)' 


To  agree  with  the  curves  of  Fig.  5 we  must  convert  the  t dependence 


where  e = to  that  of  the  figure  where  the  curves  are  plotted  for 

aVl) 

Xd  2Xd  „ ... 


a^(a=0)  a (a) 


= 2 c so 
2 2 

o f 

_ O 

S4  - "6T" 


(315) 


Kt)-5] 


In  contrast  to  the  a = 0 case  we  will  need 


(0  = 0.  2 J 1 + 


3i6eoRtHJ 


to  have  S4  =0.3  which  becomes  larger  than  the  oQ  required  for  a = 0 by  as  much 


as  a factor. 


*v/  II 

^IW1] 


for  small  e which  is  7.  25  for  e = 10  . 

o o 


Since  the  0.  9 contour  appears  to  have  an  e 7 behavior  from  the  figure, 

the  Mercier  type  extension  used  in  the  a = 0 case  will  not  work  here.  However, 

the  correspondence  of  systems-related  parameters  such  as  with  S4  breaks 

down  also.  This  is  because  the  focusing  effects  caused  by  shapes  whose  power 

-3 

spectra  decrease  more  rapidly  than  q lead  to  large  values  of  S4  from  focusing 
rather  than  from  deep  nulls. 

_2 

A check  of  the  P =10  contour  shows  it  to  follow  the  curve 
e 


<o  = 0.  9 
o 


315  ‘o*  Hi) + i. 


(3. 1.4) 


r 

i 
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(a  factor  of  4.  5 higher  than  for  the  S^  = 0. 3 curve)  closely  enough  for  our  purposes. 

_3  _2 

This  fit  becomes  very  inaccurate  for  e s 3*  10  as  the  P = 10”  contour 

-3  ® 

begins  to  approach  the  = 0. 9 contour  but  this  region  (e  s 3*  10  ) is  also  the  region 
where  the  greatest  differences  occur  between  the  a- 1 and  a/ 1 cases  so  modeling  it 
accurately  for  a=l  is  not  necessary. 

Figure  4 for  a = 0.9  illustrates  this  point.  The  S.  contours  are  similar 

-2  -3  4 

for  c s 10  and  begin  to  differ  for  e £ 3*10  . Note  the  S^  = 0. 3 contour  begins 

to  fall  off  slower  with  c approaching  the  e-dependence  found  for  the  a = 0 case.  This 

-3 

is  because  the  high  wavenumber  portion  of  the  PSD  for  a = 0.  9 has  the  q behavior  of 

-3 

the  a = 0 case.  For  « < 0.  9,  the  q behavior  appears  earlier  (lower  q)  and  one  never 
-5 

has  the  q behavior  of  a=  1 since  this  would  occur  at  a value  of  q too  low  for  the 
Bessel  functions  to  be  replaceable  by  their  asymptotic  expansion  limits. 

-3 

For  0.  9 < a<l.  0,  the  q behavior  appears  at  larger  values  of  q causing 

the  S.  = 0.  3 contour  to  differ  from  the  a=  1 behavior  at  values  of  e below  the  limit  on 

’ _2 

the  figure.  Since  both  the  a=0.  9 and  a = 1 contours  of  P =10  break  from  the 

-3  e 

expression  (3. 1. 4)  for  e 3*  10  , we  expect  the  same  type  of  break  to  occur  for 

0. 9 < « < 1. 0.  This  low  e break  will  not  be  modeled  at  this  time  for  the  0.  9 <a  < 1 
case  since  we  do  not  expect  to  be  in  this  low  e regime  when  0=1.  That  is,  values 
of  ethis  low  occur  at  early  times  after  burst  when  ft  4 1.  We  shall  therefore  approxi- 
mate the  results  of  this  section  with  the  following  functional  fit.  Let 

g(to,M1^,2^+a4^(1.aV2,fo2[,n(^)  + l] 


and  define 


®e«f  = <Po/^l  + gTFJ  ‘ 


The  limit  of  scintillation  effects  (S^  = 0. 3)  occurs  for  = 0.  2.  The  amplitude 
distribution  effectively  becomes  Rayleigh  at  = 1 with  a correlation  length 


l = 0. 6 a/<0Q  . 
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3.2  ELLIPTICAL  CONTOUR 

For  a set  of  elliptical  striations  where  the  major  axis  of  striation  i is 
at  an  angle  0.  to  the  ray  direction, 


2,  , 2#  , , 2, 2 „ Jl<qai’  „ 2,  „ 2 J2(qai) 

« («>  = T (reXV  > f (1-“  > + 2“ 

where 


a.  = a J - — 5 — - + 0 sin20 . = — — fl~+  (j52  - 1)  sin2  9. 

1 s ' -JT 


The  electron  content  of  each  striation  is  independent  of  0 and  0 . so  as  before 
a(a)  = a(0)  (1  - a2/2)-^  for  the  content  to  be  independent  of  «.  Each  striation 
then  produces  a phase  variance 


2 4,  <w2>  r . Jjji<<iai)  „ . 

% - r J aidq  ~qa~1_ 


V „ 2 J2(qai> 


) + 2a 


TqaTF 


= Iap=dJ — -2  $ 2 

01  1 + (0-1)  sin  Q. 


r Ki  = 


2 2 

The  factor  T is  the  scale  factor  for  the  q-dependence  of  tp . (q)  = <p.  (q/F  )(0  = 1). 

2 2 1 1 
That  is,  <pu  (o)  is  independent  of  0 and  0.  but  o . (q)  is  stretched  a factor  r with 

respect  to  the  0=1  case.  For  small  values  of  0^,  r >1  and  the  scintillation 

effects  will  be  worse  than  in  the  0 = 1 case. 


We  illustrate  this  by  choosing  a random  set  of  striations  for  a situation 
where  the  signal  line-of-sight  traverses  the  middle  of  a striated  region.  The 
set  of  values  of  0 . is  generated  by  choosing  values  of  x^  to  be  random  on  a line 
0 s x s64  km  while  values  of  y.  are  random  on  a line  0 sy  «s500  km.  The 
value  of  0A  is  chosen  to  be  on  the  radial  from  the  center  of  the  striations  to 
the  position  of  the  ith  striation^ 

tan9i  = i-  . 
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The  resultant  PSD  for  0 = 3 is  shown  in  Fig.  6 for  the  two  cases  a = 0, 1.  In  this 
figure  we  do  not  need  to  plot  the  loci  of  peaks  of  the  Bessel  functions  since  the 
range  of  values  of  0^  smear  out  the  far  nulls  much  as  a distribution  over  sizes 
would.  The  first  few  nulls  are  not  smeared  out  much  because  the  range 

of  0.  is  small.  The  dashed  lines  are  fits  to  the  curves  used  for  the  propaga- 

1 -3 

tion.  The  behavior  at  large  v (y=q/2ff)  will  be  proportional  to  (v/v  ) for  a = 0 

-5  ° 

and  ( v/v for  a=  1.  Since  the  PSD,  P(y),  is  normalized  to  unity  in  the  figure 

at  v = 0,  we  have  for  a = 0 


i vpVWri>l 

P ,,)  M tti[  (2irw/r)  J 


n t o c 

2ir  v a 


M 

_!  v 

M ^ 


For  a=l 


"o  = 


r.3  * 

i 

(2ir)3  ir  a 


P(y)  = 


n v a 


8J2  (2ff^a/ri) 
(2ffva/ri)2 


J_  y 

M ~ 


(w/2)/sffa 


The  curves  are  drawn  for  a (a.  =0)  = 0.  282  km  and  a (a=  1)  = 0. 282,/T  = 0.399  km 
and  show  vQ  = 0. 976  km  * and  v^  = 1. 142  km~*  so 
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ip 


= (2 ff)  = (.  976)  it  (.  282)  = 1.  595 


and 


k'A  % 

r = (ff/2)  (1.  142)  (.  399)  = 1.  567 


—J 

The  closeness  of  these  numbers  show  0 j 1 since 


_ n an/2  ,,  .2  x . 2 ,-n/2 

T.  = 0 (1  + (0  -1)  sin  0.) 


„n/2  n /„2  1N  „ 2. 

*0  (i  - ^ (/?  -i)  e.  ) , 
2 


— r 1/n 


r. 

1 


n 


0 * (1  - 0.2 ) = 1.  732  (l-4e.2) 


“T  — 2" 

which  produces  the  values  0.  = . 0198  at  a = 0 and  0.  = . 0238  at  a=l  where  the 

1 1 2 2 

difference  in  these  two  terms  are  a result  of  higher- order  terms  in  (0  - 1)  a 

which  we  dropped.  Our  procedure  for  aligning  striations  along  a radial  direc- 
tion has  led  to  the  small  value  of  8 « 8.  5*  for  the  alignment  of  major 

rms 

axes  with  the  ray  direction.  We  expect  r = 1.  6 is  sufficiently  accurate  for 
our  purposes. 

o o 

The  contours  of 


We  thus  expect  <pQ2  (0  =3)  •*  1. 6 ©o2  (0  = 1). 


S^  = 0.  3 and  0.  9 and  various  values  of  t/a  are  shown  in  Fig.  7 for  the  a = 0, 
0 = 3 case.  It  is  apparent  that  the  behavior  of  the  contours  is  similar  to 
Fig.  3 for  0=1  and  agreement  is  obtained  by  scaling  the  c parameter  by 
f 2 = 2.  56.  So  that 


S4(o0,  e) 


0 = 3 


= s4  («o,  r4.; 


0 = 1 


The  contours  of  l/a.  are  displaced  from  those  in  Fig.  3 and  have  the  form 
l/a  = 0. 6/f  (pQ)since  l must  scale  as  T1  . Comparison  of  the 
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contours  in  Fig.  8 for  the  of  1,  .8  = 3 case  with  those  shown  in  Fig.  5 indicates  that 
the  simple  scaling  agrees  sufficiently  even  though  not  as  well  as  for  the  a = 0 


case. 


The  PSD  for  8 = 10  is  shown  in  Fig.  9 for  a = 0 and  1.  From  the  curves  we 


= 1.  377  and  y.  = 1.  773  which  means 
o 1 

? % 

r.  = 2.  251  and  r.  = 2.4325. 
i l 


So  the  particular  choice  of  striations  causes  the  PSD  to  scale  as  r=  2.  25  for 


ot  = 0 and  r=  2. 43  for  a=  1.  We  have  o 


'8  = 10 


and  we  expect  the  contours 


of  S and  t/a  in  Figs.  10  and  11  to  scale  as 


s4(*o’c)  = S4 ^°o * r 

40  8 = 10  18=1 


t/a  =■ 


which  they  do. 

Because  the  scale  factors  depend  on  which  region  of  the  fireball  is 

being  scanned  by  the  line-of-sight,  we  have  chosen  a stressing  region  where 

. There  will  be  off-center  regions  where  r*  1 but  regions  with  T<1 

1 

where  for  large  8 one  has  r~  — — iTT  n will  produce  weaker  scintillation 

JJ  CSC  9i 

effects.  Regions  further  out  will  also  produce  weaker  effects  since  less  structure 
is  propagated  through  due  to  the  curvature  of  the  region  containing  the  striations. 
Even  without  the  curvature  effect,  T < 1 and  scintillation  effects  are  weaker  than 
for  circular  eross-sections  for  most  of  the  path  of  the  line-of-sight  through  the 
striations. 
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?(V) 
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3.  3 EFFECTS  FOR  A SINGLE  SIZE 

The  general  conclusions  for  scintillation  effects  are  presented  in  terms 
of 


o 


<p 

o 


eff 


Jl'+TTg 


(3.3.  1) 


where 


2 16  a % , 2 

ao  = 3E—  (reXno)  M 

(ignoring  the  slight  dependence  of  c oQ  on  a)  where  M is  the  total  number  of  striations, 
8~^  < T<  r , and 


2.2  3irfe  4n  2/9.  315  _4  2 f.  / 2 \ 3] 

g (1— o!  ) g + oi  (l'tt  /2)  g4  T c I In f ~ p J + 2 I • 

For  <peff  & 1,  the  signal  distribution  will  be  taken  as  effectively 
Rayleigh  although  the  scintillation  index,  S4,  need  not  be  unity.  In  this  case  the 
signal  autocorrelation  length 


l M 


0.6  a 


r 


(3.3.2) 


/2 


For  ^ 0.2,  S4  s0. 3 and  amplitude  scintillation  is  negligible.  The  quantity 

Xd 

2 

The  quantity  a is  defined  by  the  electron  content  of  a striation,  C s nQira  . 

2 

As  the  original  content  C = n fa  striates  and  bifurcates  the  number 


2 2 

of  striations  M = aQ  /a  so 


o o o 


= 4r  \n  a 


ra 


o ~ eA  o o * "SIT 


(3.313) 


iJAL 
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decreases  as  a decreases.  But 


s 


. 3<l-o2)  -3,  „ . a4(l-«2/2 

2 r x d + j— 2 

32a  1024  * a 


In 


/ 8?r a2  \ 

te) 


3 

+ 2 


(3.3.4) 


increases  as  a decreases. 


What  are  the  relationships  between  nQ,  a,  a,  and  r as  the  original 
striations  bifurcate?  A discussion  of  a in  Section  5 suggests  the  formula 


a(t)  = 0. 1 + \ [l9.  9 - F + /( 19.9  - F f + 0.  64  ] km 


where  t = for  t a 500  s At  t = 500  s,  a = 20  km  so  for  a = 150  km, 

o 
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M 56.  The  striations  should  be  elliptical.  j3  ~ 10  soO.  4 < T < 2.  5 depending  on 
the  region.  The  profiles  will  be  parabolic  until  n falls  to  a value  where  the  inte- 


grated conductivity  is  a few  times  the  background.  For  a 3 mho  background  and  a 

8 

30  mho  central  conductivity  at  n = 10  , the  profile  wouldn't  be  steep  until 

r 7 0 


n^  ~ 2-  10  which  corresponds  to 


= .9  . 


For  L = 10  km 


o 


<PQ(500  s)  m 108  /f  X (cm)  = 140  X (cm)  . 


For  a = 0. 9 and  d = IQ11  km 


g = 8.  5- 10'8  r2  X (cm)  + 7.  6-  10"12  r4  X 2 


T X(cm)/  1 


so 


g (r2x)~  8-io-8  r2x  r2x  < io3 


~ io~ 10  r4x2  (i  - x)  ) r2x  > io3 
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For  X = 10  cm. 


<p  = 1400  T^Jg  = 1.26  r*  r < io2 


and  the  scintillation  depends  critically  on  F (being  slight  near  the  edge  of  the 

7 

region  and  Rayleigh  in  the  rest  of  the  region).  If  n were  10  , & would  drop 

o On 

to  140  producing  negligible  scintillation  for  any  r . However  the  drop  to  10* 
means  the  profile  should  have  steepened  to  a = 0 where 


Js  T7W7i 


0eff  *7  T‘ 


with  little  benefit. 

The  fireball  region  is  thus  a region  of  high  scintillation  even  for  large 
striations  of  20  km  size  for  X = 10  cm.  The  equation  set  (3. 3. 1-3. 3.  4)  describes  the 
amplitude  scintillation  as  a function  of  the  parameters  nQ,  x,  L,  d,  a,  a , a, 
and  r (/3 ).  We  expect  a and  £ to  depend  on  nQ  , background  and  time.  We  have 
thus  far  treated  a case  where  all  striations  have  the  same  size.  We  will  investi- 
gate the  effects  of  a size  distribution  on  scintillation  in  a later  report.  The  effects 
of  a size  distribution  on  the  PSD  are  discussed  in  Section  5.  We  have  shown  in  this 
section  that  amplitude  scintillation  effects  can  be  approximated  using  analytic  fits 
to  the  PSD. 


P 
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4.  TIME -DEPENDENT  STRIA TION  DEVELOPMENT 

4. 1 DISCUSSION  OF  STRIA  TION  DEVELOPMENT  IN  HANE 

The  gradient-drift  instability  mechanism  has  been  combined  with  the 

output  of  hydrocodes  of  HANE  in  order  to  determine  those  regions  of  space  and 

time  after  burst  during  which  significant  amplification  of  irregularities  in  the 

electron  density  can  be  expected.  Generally  speaking,  these  various  studies 

have  shown  that  the  region  of  space  200  to  400  kilometers  away  from  the  center 

of  the  electron  density  distribution  tends  to  have  striation  amplifications  in  the 
2 3 

• 10-10  range  by  times  of  the  order  300  to  600  s after  burst.  Consistent  with 

; the  often  made  assumption  that  the  initial  amplitude  of  fluctuations  in  electron 

-3  -2 

\ density  are  of  order  10  - 10  times  the  unperturbed  electron  density,  the 

nominal  time  at  which  the  irregularities  in  electron  density  become  nonlinear, 
i.  e. , An  /n  of  order  1,  can  be  taken  to  be  500  seconds  after  burst  certainly 
| to  within  a factor  of  2. 

| Unfortunately,  our  present  understanding  of  the  initial  phase  of  striation 

development  is  not  sufficient  for  us  to  specify  precisely  the  form  and  scale  sizes 
that  are  likely  to  be  generated  by  the  gradient-drift  instability  mechanism.  We 
i know  that  these  scale  sizes  are  considerably  smaller  than  the  typical  50  - 100  km 

| sizes  used  in  the  hydrocode  computations.  The  scintillation  effects  to  be  expected 

j from  HANE  depend  critically  on  the  characteristics  of  the  striations.  The  pre- 

I vious  section  has  indicated  the  dependence  of  propagation  effects  on  some  of 

j these  currently  unknown  parameters.  In  the  discussion  below  we  will  describe 

those  aspects  of  striation  development  for  which  we  have  an  understanding  as 
well  as  identify  gaps  in  our  knowledge.  The  result  of  the  discussicn  provides 
the  motivation  for  the  parameterized  models  that  were  discussed  in  Section  2 and 
formed  the  basis  for  the  computations  in  Section  3. 


Figure  1(B)  provides  a schematic  illustration  of  the  electron  density 
contours  at  times  around  500  seconds  after  burst  after  the  initially  fairly  smooth 
electron  density  contours  become  highly  distorted.  The  striations  are  formed 
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FMBCSDINO  FAQS  BUkMUNOT  FILMS 


in  part  by  motion  of  regions  of  high  electron  density  to  a greater  radius  than  the 
hydrocodes  would  have  transported  them,  and  in  part,  by  motion  of  regions  of 
lower  electron  density  in  towards  the  center  of  the  electron  density  distribution. 

We  do  not  know  precisely  the  configurations  to  be  expected  at  this  time  in  the 
cylindrical  geometry  with  a radial  neutral  wind  that  is  a function  of  radial  distance. 
However,  we  summarize  below  our  current  understanding  regarding  scale  sizes, 
shapes,  and  the  time  evolution  of  striations  and  place  bounds  on  the  range  of 
values  to  be  expected. 

4. 1. 1 Initial  Development 

11  12  7 

The  analytical  and  numerical  ’ studies  of  one -dimensional  clouds 
have  some  relevance  for  the  cylindrically  symmetrical  problem.  The  unperturbed 
equilibrium  in  both  cases  has  the  Pedersen  current  flow  at  right  angles  to  the 
density  gradients  which  are  colinear  with  the  neutral  wind  velocity.  There  are, 
of  course,  geometric  factors  which  lead  to  significant  differences  on  the  global 
scale  for  the  two  cases. 

5 

The  gradient-drift  instability  results  from  disruptions  in  the  equilibrium 
current -flow-pattern  due  to  irregularities  in  the  field-line  integrated  Pedersen  (FLIP) 
conductivity,  E(r).  Thus,  in  order  to  determine  the  electric  field  pattern  that 
maintains  divergence-free  current  flow,  we  need  to  know  the  distribution  of  the 
FLIP  conductivity.  Consistent  with  the  smooth  electron  density  contours  shown 
in  Figure  1(A)  before  striation  formation,  we  may  model  the  radial  distribution 
of  the  Pedersen  conductivity  as 


2 (r)  = E Q [ 1 + x exp  (-r2/aQ2)] 


(4. 1) 


where  Eq  represents  the  FLIP  conductivity  of  the  uniform  background.  E 
for  the  ambient  ionosphere  has  typical  values  of  1 - 10  mho  depending  on  time  of 
day.  Values  of  x appropriate  for  HANE  will  be  discussed  below. 

During  the  initial  stages  when  the  fluctuation  amplitudes  are  small,  the 
development  of  the  irregularities  can  be  treated  by  linearized  equations.  The 
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problem  with  attempting  to  apply  these  results  to  HANE  is  that  the  unperturbed 
electron  density  (and  conductivity)  is  inhomogeneous  and  time-varying.  Thus, 
the  changes  in  the  electron  density  contours  are  due  in  part  to  the  description 
provided  by  the  linear ized-growth-analy sis  and  in  part  tp  changes  occurring 
in  the  evolving  electron  density  distribution.  The  latter  influences  the  validity 
of  the  linearized-growth-analysis. 

The  linearized  analysis  of  a one- dimensional  cloud  has  most  validity 

when  applied  locally.  Perkins  et  al.  ^ analyzed  the  stability  in  the  limit 

kh  » 1 where  k is  the  Fourier  wave -number  of  the  perturbed  density  in  the 

direction  of  the  current  flow  and  h = | v In  E I is  the  gradient  length  associated 

with  the  spatial  dependence  of  E (r).  The  principal  results  obtained  by  Perkins 

et  al.  of  interest  here  relate  to  the  location,  shape,  and  growth  time  of  the 

perturbed  density.  They  find  that  the  eigenfunction  maximizes  at  the  position 

r = r for  which  h is  a minimum,  h . (maximum  in  the  logarithmic  gradient), 
m m 

The  striation  is  elongated  in  the  direction  of  the  neutral  wind,  Vn,  with  a length 
that  is  8 = (2khm)^  times  the  transverse  width,  l/k.  The  time  scale  for  r-monen- 
tial  growth  is 


T 

g 


<hm/Vs)[l+<lAhm» 


(4.2) 


where  Vg  is  the  magnitude  of  the  slip  velocity  of  the  neutrals,  Vg  = VR  - 
transverse  to  the  magnetic  field.  Because  diffusion  coefficients  are  generally 
small,  diffusion  is  not  likely  to  be  a dominant  process  until  late  times  when 
large  density  gradients  have  formed. 

Table  4. 1 shows  some  of  the  characteristic  numbers  that  result  when 
the  conclusions  from  the  linearized  analysis  are  applied  to  the  cases  of  both 
planar  clouds,  shown  above  the  dashed  line,  and  cylindrical  clouds  associated 
with  HANE,  shown  below  the  dashed  line.  The  cases  above  the  dashed  line  are 
chosen  to  represent  the  parameters  associated  with  the  numerical  simulations 
of  NRL6’  7 when  applied  to  clouds  having  a conductivity  profile  given  by  Eq.  (4. 1). 
The  third  column  gives  the  position  of  the  maximum  in  the  eigenfunction  pre- 
dicted by  linear  theory.  We  see  that  for  higher  conductivity  ratios  this 
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Table  4. 1 Cloud  characteristics  associated  with  initial  striation  development 
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maximum  relative  deformation  occurs  further  out  from  the  region  of  peak  density 
toward  the  tail  of  the  Gaussian.  Column  4 gives  the  value  of  the  cloud  conductivity 
in  units  of  at  the  position  of  the  maximum.  Even  for  clouds  with  very  high  peak 
conductivity  ratios,  the  density  contours  that  deform  most  rapidly  correspond  to 
a value  of  cloud  conductivity  of  order  the  background  conductivity. 

Column  5 lists  the  minimum  value  of  the  logarithmic  density  gradient 
length  that  occurs  for  each  of  the  cases.  Column  6 indicates  the  characteristic 
full-width  of  the  striations.  In  the  planar  cloud  cases  above  the  dashed  line,  the 
values  given  are  obtained  from  the  observation  that  in  the  12  kilometer  dimension 
shown,  there  were  5 and  3 striations  growing  initially,  in  the  two  cases  respectively. 
This  dimension  is  also  consistent  with  the  full-width  after  the  striations  have 
become  nonlinear.  In  the  HANE  cylindrical  cloud  case  shown  below  the  dashed 
line  we  have  just  given  characteristic  numbers  associated  with  a striation  of  a 
characteristic  radius  a = 20  kilometers.  The  seventh  column  shows  that  these 
numbers  imply  that  khm  » 1 . The  eighth  column  provides  an  estimate  of  the 
number  of  such  20  kilometer  size  striations  that  might  exist  around  the  cylindri- 
cal radius  of  r = r . The  characteristics  growth  time,  r ~ h /V  is  shown  in 
column  9 where  for  the  planar  cloud  case  a velocity  of  0.  1 km/s  was  used  while 
in  the  HANE  case  a typical  velocity  of  1 km/s  was  used. 

We  now  compare  the  numbers  relevant  for  the  linearized  analysis  shown 
in  Table  4. 1 with  numerical  simulations  of  the  planar  cloud  case.  Figures  12 
(from  Ref.  7)  and  13  (from  Ref.  6)  illustrate  the  x - 1 (at  t = 480  s)  and  x = 10 
(at  t = 240  s)  cases,  respectively.  In  each  of  these  figures  the  maximum  distor- 
tions of  the  equal  density  contours  are  comparable  to  the  size  of  the  striations 
which  is  a criterion  that  the  distortions  are  ending  the  linear  phase  and  entering 
the  nonlinear  phase.  Note  that  the  positions  of  the  maximum  distortions  are 
located  at  a distance  further  from  the  maximum  density  in  the  case  of  the 
higher  conductivity  ratio.  At  the  times  indicated  in  column  10  of  Table  4.  1, 
the  initial  fluctuation  amplitudes  of  5%  and  3%,  respectively,  will  have  been 
amplified  by  the  factor  A shown  in  column  11.  This  amplification  is  sufficient 
that  the  linearized  theory  is  no  longer  applicable  and  the  striation  development 
has  entered  the  nonlinear  regime  consistent  with  the  distortions  shown  in  the 
respective  figures. 
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Figure  12.  Iso-Pedersen  conductivity  contours  of  a one -dimensional 

plasma  cloud  with  x = 1 at  t = 480  s after  initial  perturbation 
amplitudes  of  ~ 5%  (taken  from  Ref.  7). 
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Figure  13.  Iso-Pedersen  conductivity  contours  of  a one -dimensional 

plasma  cloud  with  x = 10  at  t = 240  s after  initial  perturbation 
amplitudes  of  ~ 3%  (taken  from  Ref.  6:  R = ® case). 
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At  late  time  in  the  nonlinear  regime,  the  above  description  is  not  appro- 
priate. In  that  regime  the  striations  do  not  increase  in  amplitude  exponentially 
with  time  but  maintain  a nearly  constant  electron  density.  The  striations  move  with 
a velocity  that  is  dependent  on  the  local  slip  velocity  and  the  ratio  of  the  FLIP 
conductivity  of  the  striation  to  the  FLIP  conductivity  in  the  background  region  in  the 
vicinity  of  the  striation.  Furthermore,  these  isolated  striations  undergo  additional 
deformations  that  are  sometimes  simply  described  as  bifurcations.  We  will  discuss 
the  time-scale  for  these  bifurcations  in  Section  4. 1. 2 below. 

We  do  not  have  an  adequate  understanding  of  the  transition  between  the 
linear  regime  and  the  nonlinear  regime  except  in  highly  idealized  cases.  Crucial 
questions  regarding  the  scale  size  of  the  striations  and  the  time  scale  for  their 
development  require  further  study. 

If  an  initially  planar  cloud  is  seeded  with  perturbation s having  a periodicity 

11  12 

equal  to  2 tr/k,  then  both  linear  and  nonlinear  analyses  show  that  the  resulting 

striations  maintain  this  periodicity.  In  the  linearized  regime,  the  perturbation 

density  grows  exponentially  in  time  with  a time  constant,  r , close  to  that  given 

6 

by  Eq.  (4.  2).  However,  linear  theory  so  far  has  been  unable  to  specify  the  eventual 
scale  size  of  striations  in  the  nonlinear  regime.  A naive  interpretation  of  the 
linear  theory  would  indicate  that,  according  to  Eq.  (4.2),  the  smallest  striation 
scale  size  would  grow  the  most  rapidly.  It  would  also  indicate  that  for  the  case 
of  large  conductivity  clouds  the  striation  growth  would  affect  primarily  the  low 
density  regions  of  the  cloud,  near  r = rm  where  rm  is  the  location  of  the  maxi- 
mum in  the  eigenfunction. 

However,  as  the  striations  approach  the  nonlinear  regime  neither  of  the 
above  conclusions  is  correct.  The  small  scale  perturbations  which,  according  to 
linear  theory,  grow  the  fastest  also  "saturate”  the  fastest.  Larger  scale  striations 
which  continue  growing  then  dominate  the  smaller  scale  structures.  The  influence 
that  the  value  of  the  cloud  conductivity  ratio,  x , the  density  profile  of  the  cloud, 
connections  to  other  layers  of  ionization  on  the  same  magnetic  field  line,  and 
possibly  diffusion,  have  on  the  important  problem  of  specifying  the  scale  size  of 
nonlinear  striations  has  not  yet  been  determined. 
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In  regard  to  the  striation  scale  sizes  that  we  might  expect  from  HANE 

from  the  available  examples,  it  would  appear  that  it  is  reasonable  that  the  initial 

values  would  be  in  the  range  3 < kh  < 30.  This  result  would  give  initial  char- 

m 

acteristic  striation  half-widths,  a - ff/2k,  in  the  range  from  several  kilometers 
to  several  tens  of  kilometers.  Furthermore,  an  examination  of  column  11  of 
Table  4. 1 indicates  that  by  times  of  around  500  s following  burst,  the  amplification 
factors  are  generally  in  the  range  consistent  with  code  calculations.  Both  the 
simplified  model  and  the  more  extensive  code  calculations  indicate  that  the  initi- 
ally linear  striations  have  become  nonlinear  by  approximately  500  s. 

4. 1.  2 Nonlinear  Development 

We  now  turn  to  a discussion  of  some  of  the  expected  characteristics  of 
striations  when  they  are  in  the  nonlinear  phase.  For  the  case  of  HANE  we  assume 
that  the  striations  have  reached  the  nonlinear  phase  at  time  s greater  than  500  s 
after  burst.  We  will  discuss  some  concepts  related  to  the  initial  sizes  of  the 
striations,  the  shapes  of  the  striations,  the  density  profiles  within  the  striations, 
and  the  time  scale  for  their  development.  A discussion  of  the  time  development 
of  striations  is  given  in  the  next  section. 

Once  the  ionization  in  HANE  has  reached  the  stage  shown  schematically 
in  Figure  1(B),  we  can  apply  our  understanding  of  the  processes  that  affect 
the  evolution  of  these  striations  which  we  will  characterize  as  having  a 
typical  scale  size  a.  In  discussing  these  concepts  we  will  not  include  the  mutual 
interactions  between  neighboring  striations.  This  neglect  will  modify  some  of 
the  details  but  will  not  affect  the  general  concepts.  We  have  already  indicated 
that  the  probable  size  of  these  striations  is  determined  by  processes  not  fully 
understood  during  their  initial  stages  of  development.  Indications  are  that  for 
the  case  of  HANE  they  are  probably  in  the  range  of  several  kilometers  to  several 
tens  of  kilometers  in  scale  size. 

The  conductivity  associated  with  a HANE  cloud  is  an  important  para- 
meter for  describing  its  development.  During  the  first  tens  of  seconds  following 
a high  altitude  explosion  the  concept  of  conductivity  is  inappropriate  because  the 
magnetic  field  lines  are  distorted  and  the  electrostatic  approximation  is  not  valid. 
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If  the  concept  were  applied  unthinkingly,  however,  the  evaluation  of  the  FLIP  con- 
ductivity would  probably  produce  values  in  the  range  of  many  of  hundreds  to  thousands 
of  mho.  However,  this  value  would  rapidly  decrease  with  time  because  the  neutral 
density  decreases  as  a result  of  being  heaved  to  higher  altitudes. 

The  contribution  to  the  FLIP  conductivity  obtained  by  integrating  over  a 
small  height  range  A , is  given  by 


« ^ - eny  A 

5 L ~m  ire 


(4.3) 


where  e,  n,  B,  u,  ft,  and  0 are,  respectively,  the  electronic  charge,  ion  number 
density  (equal  to  the  electron  concentration),  magnetic  field  strength,  ion-neutral 

collision  frequency,  ion  gyrofrequency,  and  magnetic  field  dip  angle.  By  using 

-1  13 

B = 0.  5 gauss,  O = 300  s appropriate  for  ions  of  mass  16,  and  v = 2 x 10  p 

3 

where  p is  the  neutral  mass  density  expressed  in  g/cm  , (consistent  with  an 
approximate  effective  collision  cross  section  of  ~ 5 x 10~15  cm\  and  0~  70? 


we  obtain 


6 £ ~ 2. 3 x 10  nPA  mho 


(4.4) 


where  n and  A are  expressed  in  units  of  cm  and  km,  respectively.  As  an 
example,  the  FLIP  conductivity  of  a plasma  cloud  produced  by  a 48  kg  barium 

7 _3 

release  at  ~ 10  minutes  after  release  when  n ~10  cm  , A ~ 25  km,  and 
-13  3 

p ~ 4. 5 x 10  g/cm  , is  ~ 26  mho  when  Eq.  (4.  4)  is  used.  This  value  is  in 
excellent  agreement  with  more  detailed  estimates  given  in  Ref.  13. 

Accurate  values  of  the  FLIP  conductivity  of  a HANE  cloud  as  a function 
of  time  can  be  obtained  from  late-time  electrostatic  codes  such  as  MELT,  but 
these  values  are  not  currently  part  of  the  output.  It  is  difficult  to  make  accurate 
estimates  by  using  Eq.  (4.  4)  because  both  n and  p vary  by  many  orders  of 
magnitude  along  a magnetic  field  line  and  this  variation  changes  with  time.  How- 
ever, we  have  attempted  to  use  the  expression  given  in  Eq.  (4.  4)  to  get  crude 
estimates  of  the  range  of  expected  FLIP  conductivities.  We  have  examined 
typical  outputs  of  the  MELT  code  at  times  of  500,  1000,  1500,  and  2000  seconds 
after  a burst.  Because  many  neutrals  have  been  heaved  to  high  altitude,  it 
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appears  that  at  500  s,  the  FLIP  conductivity  ~ 100  mho.  By  1000  s,  it 
appears  that  this  value  may  have  decreased  by  an  order  of  magnitude  or  more. 

For  the  time-scale  500  - 2000  s after  burst  the  FLIP  conductivity  of  the 
large  cloud  of  ionization  created  by  a high  altitude  nuclear  burst  is  of  order 
10  - 100  mho. 

Forty-eight  kilogram  barium  releases  at  altitudes  of  180  to  190  km 

13 

produce  clouds  of  ionization  that  have  a FLIP  conductivity  of  order  30  mho  , 
and  a typical  scale  size  a ~ 3 km.  Thus,  it  is  reasonable  to  expect  that  some 
of  the  observed  behavior  of  barium  releases  conducted  in  the  ionosphere  at  dusk 
may  provide  appropriate  descriptions  of  the  subsequent  development  of  the  large 
scale  striations  shown  in  Figure  1(B),  particularly  at  later  times.  The  scale 
size  of  a 20  km  striation  in  HANE  is  approximately  a factor  of  7 larger  than  a 
3 km  barium  cloud.  However,  at  late  times  in  the  1000  s time  range  following 
a burst,  the  neutral  wind  velocity  at  several  hundred  kilometers  radius  from  the 
burst  field  line  is  of  order  500  m/s  which  is  a factor  of  7 larger  than  a typical 
ionospheric  neutral  wind  velocity  of  70  m/s.  Thus,  the  time  scale  for  the  sub- 
sequent development  of  the  large-scale  nonlinear  striations  that  can  be  assumed 
to  be  present  at  500  seconds  after  burst  is  in  the  same  range  as  the  time  scale 
for  the  development  of  smaller-scale  ionospheric  barium  clouds.  Each  of  the 
clumps  shown  in  Figure  1(B)  is  likely  to  develop  in  a manner  very  similar  to 
typical  barium  clouds  such  as  Spruce  or  those  released  during  the  STRESS  test 
series. 

Barium  releases  conducted  in  Florida  during  winter  behaved  as 
though  the  ratio  of  cloud  to  background  FLIP  conductivity  were  in  the  range  of 
3 to  6,  so  we  estimate  the  effective  ionospheric  conductivity  to  be  5 to  10  mho. 
These  values  are  significantly  greater  than  is  expected  from  a nighttime  iono- 
sphere which  should  be  in  the  range  of  1 mho.  The  larger  value  is  explained 
by  the  hypothesis  that  the  polarization  electric  fields  induced  by  the  presence 
of  the  barium  cloud  is  projected  along  the  magnetic  field  lines  into  the  daytime 

(in  the  case  of  Spruce)  or  late  afternoon  (in  the  case  of  the  STRESS  test  series) 

14 

ionosphere  in  the  southern  conjugate  region  . Thus,  the  appropriate  background 
ionospheric  conductivity,  S , to  use  in  Eq.  (4. 1)  may  range  from  a low  of  the 
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order  of  1 mho  up  to  a high  in  the  10-20  mho  range  corresponding  to  the  presence 
of  nighttime  or  daytime  ionospheres.  The  appropriate  value  of  x in  Eq.  (4. 1)  for 
HANE  may  range  from  < 1 at  late-time  for  a daytime  burst  to  possibly  as  large  as 
> 100  at  500  s for  a nighttime  burst.  The  strong  suggestion  is  made  that  rather  than 
providing  just  a simulation  of  gradient -drift  physics,  observations  of  barium  clouds 
may  provide  a direct  visualization  of  the  equivalent  development  of  structure  for 
some  cases  in  HANE. 

We  now  discuss  various  characteristics  of  the  striations  in  the  nonlinear 
regime.  The  first  characteristic  is  the  peak  electron  concentration  in  the  striations. 
Striations  are  formed  by  the  tendency  for  the  regions  of  large  conductivity  to  become 
more  coupled  to  the  neutral  wind.  Thus,  the  peak  electron  concentrations  in  stri- 
ations are  essentially  comparable  to  the  peak  electron  concentrations  that  are  pro- 
vided by  the  output  of  the  large-scale  codes.  Many  of  the  striations  will  have 

8-3 

electron  concentrations  of  order  10  cm  for  thousands  of  seconds. 

The  shape  and  profile  of  the  striations  depend  on  the  conductivity  ratio 
x.  We  have  already  discussed  the  uncertainty  in  our  knowledge  of  the  initial  ale 
size  of  the  striations,  but  as  we  discussed  in  the  previous  subsection,  we  expect 
values  of  khm  in  the  range  of  1 - 100.  During  the  first  500  s when  we  expect  the 
conductivity  of  the  burst-produced  ionization  to  be  very  large,  an  examination  of 
the  tabulations  given  in  Table  4. 1 suggests  initial  scale  sizes  of  striations  from 
kilometers  to  several  tens  of  kilometers.  Associated  with  this  range  of  striation 
scale  sizes,  we  expect  an  initial  ratio  of  length  to  width,  defined  earlier  as  the 
parameter  /3,  in  the  range  from  1 to  10.  Two  sets  of  arguments  lead  to  this 
conclusion.  First,  the  linear  theory  provides  an  estimate  /9  ~ (2  khm)^  when 
khm  » 1 which  is  consistent  with  1 < 0 < 10.  So,  from  linear  theory  one  expects 
elongated  initial  striations.  The  second  argument  follows  from  an  analysis  of  the 
polarization  electric  field  pattern  associated  with  nonlinear  striations.  For 
large  conductivity  ratios,  the  rate  of  steepening  of  the  original  high  density  con- 
tours in  the  direction  of  the  neutral  wind  is  much  lower  than  for  moderate  con- 
ductivity ratios.  Thus,  after  the  nonlinear  striations  have  been  formed,  the 
typical  length  in  the  direction  of  the  neutral  wind  may  be  a reflection  of  the 
original  scale  size  in  that  direction,  i.  e. , a dimension  of  order  150  km.  If.  indeed. 
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the  dimension  of  the  striations  in  the  radial  direction  were  this  large,  then  it 
is  easy  to  imagine  values  of  j3  ~ 10. 

Now  we  discuss  the  profile  of  the  electron  density  within  a striation. 

A principal  feature  of  plasma  cloud  dynamics  is  the  steepening  of  the  backside 
of  a plasma  cloud.  However,  this  steepening  phenomenon  is  most  prominent 
for  those  parts  of  the  cloud  or  striation  that  have  values  of  conductivity  less 
than  several  times  the  background  conductivity  2 . Conclusions  drawn  from 
both  linear  and  nonlinear  theory  support  this  claim.  From  linear  theory  we 
pointed  out  when  we  discussed  column  4 in  Table  4. 1 that  the  values  of  the 
conductivity  at  which  the  linearized  perturbation  amplitude  is  largest  occurs 
for  those  values  of  conductivity  that  are  only  a few  times  the  background  value 
even  when  the  peak  conductivity  is  extremely  large,  i.  e. , relatively  low-valued 
density  contours  undergo  initial  maximum  deformation. 

Arguments  from  nonlinear  analyses  rely  on  the  realization 
that  steepening  occurs  as  a result  of  shear  flow.  In  the  electrostatic  approxi- 
mation that  is  applicable  for  the  description  of  striation  development,  the  shear 
in  the  flow  exists  because  of  gradients  in  the  electric  field.  Gradients  in  the 
electric  field  result  from  the  presence  of  polarization  charges  whose  magnitude 
and  location  are  determined  by  Poisson's  equation.  In  the  simplified  one -layer 
model,  the  requirement  that 

v • J = 0 where  J^SfE  + V^xB)  (4.  5) 

is  the  field-line- integrated  current  density  across  B can  be  written  as 

v • E = - (E  + Vn  x §)  • v In  L (r)  . (4.6) 

This  equation  indicates  that  the  extremes  in  the  gradients  of  E are  at  the  posi- 
tions where  h = | v In  2 | * is  a minimum.  An  examination  of  column  4 of 
Table  4. 1 shows  that  even  for  large  x , this  position  occurs  at  relatively  low 
values  of  conductivity.  Hence,  we  expect  those  portions  of  the  striations  that  have 
FLIP  conductivities  only  a few  times  the  background  to  develop  large  density 
gradients.  However,  if  the  conductivity  ratio  associated  with  the  peak  electron 
concentration  is  of  order  10  or  more,  we  expect  the  steepening  rate  of  the  high 
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density  portions  of  the  striation  to  be  considerably  less  than  that  of  the  low  density 

,.  15 

portions  . 
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Numerical  solutions  of  the  nonlinear  equations  ’ , corresponding  to 
the  cases  x = 1 and  10,  support  these  comments.  For  the  low  conductivity  case 
considered  in  Ref.  7,  the  density  gradients  developed  at  late  times  are  signifi- 
cantly larger  than  the  original  density  gradients  present  initially.  However, 
the  late -time  density  contours  in  Ref.  6 for  Re  = ® indicate  that  large  density 
gradients  are  developed  only  on  the  outermost  density  contours.  The  density 
gradients  in  the  direction  of  the  neutral  wind  in  the  vicinity  of  the  peak  conduc- 
tivity are  similar  to  those  present  initially,  and  have  not  steepened  appreciably. 

Hence,  in  the  case  of  striations  whose  conductivity  ratios  exceed  5t  it  is  appro- 
priate to  describe  the  density  profile  as  having  a small  gradient  until  the  conduc- 
tivity falls  to  several  times  the  background  conductivity  at  which  point  the  density 
gradients  become  exceedingly  large.  These  remarks  refer  primarily  to  gradients 
in  the  direction  of  the  neutral  wind.  In  both  high  and  low  conductivity  cases,  however, 
very  large  density  gradients  are  developed  in  the  direction  transverse  to  the  direction 
of  elongation  because  this  is  the  region  where  the  shear  flow  is  the  greatest. 


One  might  suspect  that  diffusion  would  prohibit  the  development  of  very 
large  density  gradients.  However,  a simple  plasma  model  calculation  shows  that  this 
conclusion  is  not  correct.  In  fact  we  shall  show  that  very  large  gradients  in  the  low 
density  region  of  a striation  can  still  exist. 


The  continuity  equation  for  incompressible  flow  with  diffusion  added  can 

be  written 


dn  _ 

+ v-vn  = v-Dvn  (4.7) 

St  e e 


where 


K(T  + T.)  v 


is  the  classical  ambipolar  diffusion  coefficient  for  diffusion  perpendicular  to  the 
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magnetic  field,  w = 9 x 10  s"  is  the  electron  cyclotron  frequency  and 
ce 

v - v + v ■ is  the  sum  of  the  electron  collision  frequencies  for  collisions 
e en  ei 

with  neutrals  and  ions,  respectively.  With 


v ~ 2 x 10  ® n 
en  o 


-3 

v ■ ~ 1 x 10  n 
ei  e 


(4.9) 


appropriate  for  T ~ T.  ~ 0. 1 ev,  we  find  that  v.>v  for  n > 2 x 10  n . 
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For  neutral  particle  densities  less  than  4 x 10  cm-3  corresponding  to 
-13  3 

p ~ 10*  g/cm  , electron-ion  collisions  dominate  electron-neutral  collisions 
5-3 

whenever  n > 10  cm  . In  this  case  the  diffusion  coefficient  is  proportional 
e 


to  ng  and  it  can  be  written  approximately 


D = Dn  ~ 4 x 10  ^ n m^/s 
e e 


(4.  10) 


with  n in  units  of  cm  . 
e 

We  wish  to  examine  the  density  profile  at  the  tip  of  the  striation  and 
consider  density  gradients  in  the  x-direction  alone.  In  the  frame  moving  with  the 
tip  of  the  striation,  the  flow  in  the  low-density  ionization  exterior  to  a rod-like 
striation  has  a stagnation  point  at  which  we  take  the  origin  of  our  coordinate 
system.  In  this  frame,  the  incompressible  exterior  flow  field  can  be  approximated 
near  the  stagnation  point  by 


v =-coxe  + uye 

w ^ y 


(4.  11) 


in  this  planar  model.  The  velocity  shear  co  has  a value  of  order  V /a.  By 

s 

combining  Eqs.  (4. 10)  and  (4. 11)  with  Eq.  (4.  7),  we  look  for  a profile  of  the 
electron  density  that  satisfies  the  steady-state  continuity  equation  (4.  7), 


D -r~ 
dx 


+ 0>x 


= 0. 


(4.  12) 


We  have  not  found  an  analytic  solution  of  this  nonlinear  inhomogeneous  equation, 

but  we  have  found  that  there  is  a solution  for  which  n -0  as  x-x  that  can 

e o 


be  developed  as  a power  series  expansion; 


(4.  13) 


where  A is  an  arbitrary  constant.  This  solution  shows  that  v ng  -*  « as 

x - xq  and  ng  - 0.  The  constant  A and  the  value  of  xq  are  determined  by  boundary 

conditions  whose  application  requires  a knowledge  of  the  full-solution  and  depends  on 

the  value  of  n at  x = 0 and  the  striation  scale  size  a.  The  scale  size  for  the 
e 

variation  in  ng  at  low  values  of  density  is 


x 


o 


3 m 


(4.  14) 
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for  electron  concentrations  ~ 10  cm  and  striations  with  a - 1 km  and  V ~ 500  m/s. 

s 

This  small  value  for  scale  length  compared  to  a justifies  the  use  of  an  abrupt  fall-off 
in  density  as  characterized  by  the  parameter  a in  Section  2. 


In  summary,  the  discussions  in  this  section  have  provided  the  physical 
basis  for  the  models  of  striation  shapes  and  density  profiles  that  were  parameterized 
in  Section  2.  It  would  appear  that  values  of  a and  p in  the  ranges  Osasl  and 
1 s p «£  10  span  the  values  that  are  appropriate  for  various  bursts  at  certain  times 
during  the  evolution  of  the  striations.  The  largest  uncertainty,  in  addition  to  those 
represented  by  these  ranges  of  values,  has  to  do  with  the  scale  size  of  striations  to 
be  expected.  In  the  next  section,  we  will  discuss  a model  for  the  time  dependence 
of  the  changes  in  the  sizes  and  numbers  of  striations  after  the  first  500  s following 
a nuclear  burst. 


4. 2 DISCUSSION  OF  TIME-DEPENDENT  STRIATION  DEVELOPMENT 


It  is  well  known  that,  when  a plasma  cloud  is  placed  into  the  ionosphere, 
it  eventually  develops  structure  that  at  late  times  has  the  appearance  erf  isolated 
rods,  called  striations,  aligned  with  the  earth's  magnetic  field.  A one -dimensional 
cloud  that  is  elongated  in  the  direction  of  the  current  flow,  so  that  3 • 7E=  0,  is 
in  equilibrium.  If  there  were  no  fluctuations  in  the  electron  density  it  would  not 


break  up  into  smaller-scale  striations.  However,  this  equilibrium  configuration 
is  unstable.  If  the  cloud  were  seeded  with  small  amplitude  perturbations  in 
electron  density,  these  perturbations  would  initially  grow  in  amplitude  on  the 
back  side  (the  side  out  of  which  the  neutral  particles  are  flowing).  If  the  initial 
seed  enhancements  were  periodic,  the  developing  structure  would  have  the 
appearance  of  a periodic  distribution  of  sheets  (see  Fig.  14  taken  from  Ref.  12). 

If,  on  the  other  hand,  the  perturbations  in  electron  density  were  randomly  placed, 
or  contained  a large  range  of  scale  sizes,  then  sheet-like  protuberances  of 
different  size  and  amplitude  would  begin  to  grow  (see  Fig.  12  from  Ref.  7). 

The  striations  form  as  a result  of  the  tendency  for  the  regions  of  maximum  con- 
ductivity to  move  with  a velocity  closer  to  that  of  the  neutral  wind. 

As  the  striation  continues  to  develop  it  continues  to  have  a sheet-like 
appearance.  However,  the  maximum  in  the  conductivity  tends  to  approach  the 
tip  region  of  the  striation  and  the  ionization  that  is  left  behind  has  considerably 
lower  values  of  conductivity.  Since  the  ionization  flow  in  the  plane  transverse  to 
the  magnetic  field  is  incompressible,  the  width  of  the  lower  conductivity  region 
of  the  sheet  becomes  quite  thin.  This  fact  explains  the  rod-like  appearance  of 
the  region  of  higher  electron  density  concentrated  near  the  tip  when  viewed  across 
magnetic  field  lines. 

There  is  no  equilibrium  configuration  for  a tw)  -dimensional  cloud  with 
finite  conductivity  gradients.  The  first  manifestation  of  the  distortion  of  the 
electron  density  contours  is  a steepening  of  the  backside  as  the  maximum  in  con- 
ductivity moves  with  a velocity  closer  to  that  of  the  neutral  particles  compared 
to  the  velocity  of  the  lower  conductivity  regions.  An  example  of  these  distortions 
is  shown  in  Fig.  15  from  Ref.  12.  If  two  perturbations  in  electron  density  are 
placed  on  an  otherwise  cylindrically  symmetrical  two-dimensional  cloud,  each 
maximum  will  tend  to  form  its  own  sheet  and  the  cloud  initially  breaks  into  two 
pieces.  This  process  is  called  bifurcation  and  the  classical  shapes  that  result 
are  shown  in  Fig.  16  taken  from  Ref.  12.  This  bifurcation,  or  budding,  process 
appears  to  be  a common  stage  in  the  evolution  of  the  tips  of  striations.  These 
classical  shapes  are  often  seen  to  develop  in  the  photographs  of  barium  clouds 
taken  up  the  field  lines  and  in  numerical  simulations  of  late-time  striations. 


Isodensity  contours  and  striations  for  a one-dimensional  cloud  with  a periodic 
initial  perturbation  (taken  from  Ref.  12). 


Figure  15.  Isodensity  contours  for  an  initially  circular  plasma  cloud  with  no 
initial  perturbation  showing  backside  steepening  (taken  from  Ref. 
12). 


The  process  by  which  an  initially  cylindrically- symmetric  plasma  cloud 
breaks  up  into  smaller-scale  striations  is  not  certain.  One  viewpoint  suggests  that 
after  backside  steepening  has  destroyed  cylindrical  symmetry,  the  cloud  bifurcates, 
producing  smaller-scale-size  structures  which  then  continue  to  evolve  by  a series  of 
bifurcation  processes.  Another  view  tends  to  suggest  that  the  backside  steepening 
process  produces  a more  one-dimensional-like  cloud  on  the  backside  which  becomes 
more  unstable  than  an  initially  cylindrical  cloud.  This  viewpoint  would  lead  to  the 
conclusion  that  this  more  one-dimensional -like  cloud  would  break  up  into  a series 
of  sheets  whose  scale  lengths  are  determined  by  processes  that  have  not  yet  been 
clearly  identified.  There  is  no  doubt  that  these  sheets  continue  to  evolve  by  a 
bifurcation  (or  budding  or  cleavage)  process. 

No  numerical  simulation  of  an  initially  cylindrically-symmetric  cloud  in 
a one-layer  ionospheric  model  has  spontaneously  generated  smaller  scale  structures 
by  either  developing  a series  of  sheets  or  a bifurcation  process.  It  is  not  known  if 
the  one- layer  model  would  cause  a cylindrical  cloud  to  eventually  bifurcate 
or  striate.  A recent  study  by  Ossakow  et  al.  L leads  to  the  suggestion  that  the  back- 
side steepening  of  a cloud  is  perhaps  the  first  stage  in  the  beginning  of  a bifurcation. 
However,  in  the  process  of  evolving, the  steep  backside  becomes  quasi-one-dimensionai. 
It  would  appear  that  a number  of  sheets  could  form  in  this  region,  and  their  develop- 
ment could  result  in  additional  sheets  being  formed  of  smaller  size.  At  late  times, 
each  of  these  sheets  will  have  elongated  in  the  direction  of  the  neutral  wind  and  will 
have  the  appearance  of  sheets  with  rod-like  tips.  This  sheet-like  appearance  is 
similar  to  that  seen  in  photographs  of  barium  clouds  taken  at  late  times  up  the 
magnetic  field  line.  Specific  examples  are  the  Spruce  and  Redwood  ion  clouds 
released  during  test  series  SECEDE  and  the  Carolyn  ion  cloud  released  during  the 
STRESS  test  series. 
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A numerical  study  by  Doles  et  al.  suggests  that  interactions  with  other 
layers  in  the  ionosphere  can  lead  to  the  spontaneous  formation  of  sheets  without  the 
need  for  initial  large  seed  amplitudes.  At  the  moment  it  is  difficult  to  say  what  the 
dominant  processes  are  that  lead  to  the  generation  of  nonlinear  sheets  in  large 
plasma  clouds  in  the  ionosphere.  Each  of  the  mechanisms  mentioned  above  may 
have  some  role  to  play.  The  identification  of  the  role  that  each  of  those  factors 
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plays  will  be  necessary  before  one  can  specify  with  certainty  the  mechanism  that 
determines  the  scale  length  of  the  nonlinear  striations. 

We  have  already  alluded  to  the  fact  that  the  large  scale  striations, 
likely  to  exist  at  around  500  s following  HANE,  may  have  an  appearance  somewhat 
similar  to  those  shown  in  the  schematic  diagram  in  Fig.  1(B).  We  expect  each 
of  these  large-scale  striations  to  continue  to  evolve  in  a manner  very  similar  to 
ionospheric  plasma  clouds  for  two  reasons.  First,  the  typical  ratio  of  the  FLIP 
conductivity  of  the  cloud  to  that  of  the  background  ionosphere  is  expected  to  be  in 
the  range  from  < 1 to  100  which  overlaps  the  range  of  values  of  x associated  with 
many  ionospheric  barium  releases.  Second,  by  500  s the  regions  of  maximum 
conductivity  are  probably  located  in  the  range  of  from  200  to  400  kilometers  away 
from  the  center  of  the  burst  so  that  the  radial  neutral  wind  tends  to  be  flowing 
through  each  of  the  clumps  of  ionization  in  one  direction  similar  to  the  wind  flow 
through  an  ionospheric  barium  plasma  cloud. 

The  equations  (without  diffusion)  that  describe  the  behavior  of  iono- 
spheric plasma  clouds  do  not  contain  any  fundamental  scale  lengths,  but  do  have 

the  neutral  slip  velocity,  V , as  a fundamental  velocity.  The  only  lengths  that 

s 

enter  the  problem  are  those  associated  with  the  initial  distribution  of  conductivity. 
According  to  Table  4. 1,  two  such  lengths  associated  with  Gaussian  distributions 
are  the  Gaussian  radius  aQ  and  the  maximum  gradient  length  h^.  It  is  not 
known  what  effect  the  assumed  initial  density  profile  has  on  the  future  develop- 
ment of  such  a cloud.  However,  all  clouds  with  the  same  shape,  density  profile, 
and  peak  conductivity  should  evolve  on  a time  scale  that  is  proportional  to  a 
characteristic  size,  say,  ao«  Thus,  after  striations  have  become  fully -developed 
and  the  maximum  in  the  conductivity  has  concentrated  near  the  tip  so  it  is  some- 
what rod-like,  the  tip  of  the  striation  should  continue  to  evolve  in  a manner 
precisely  similar  to  the  manner  in  which  the  original  cloud  evolved,  except  on 
a time  scale  somewhat  faster  because  of  the  smaller  size.  This  concept  leads 
to  a rapid  cascade  producing  larger  numbers  of  striations  of  smaller  and  smaller 
size.  Before  continuing  the  discussion  of  the  effect  of  size  on  the  development 
time  of  striations,  we  first  indicate  the  dependence  of  the  development  time  on 
the  shape  of  the  cloud  and  the  conductivity  ratio  x . 
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Even  though  we  cannot  specify  precisely  the  method  by  which  a large 
plasma  cloud  breaks  up  into  smaller-scale  striations,  we  can  indicate  the  scaling 
of  the  time  dependence  on  the  shape  of  the  cloud  and  its  conductivity.  In  a frame 
moving  with  an  appropriate  field-line  averaged  neutral  wind  velocity,  the  zeroth- 
order  equation  for  current  continuity  has  the  simple  form 


v • 2(r)  E ’ = 0 


(4. 15) 


where  E ' is  an  effective  electric  field  . A plasma  cloud  represents  an  enhance- 
ment of  the  conductivity  in  a region  where  the  neutral  wind  is  blowing  through  back- 
ground ionization.  The  gradients  in  the  conductivity  generate  polarization  electric 
fields  that  allow  the  cloud  to  move  across  the  magnetic  field  lines  with  part  of  the 
velocity  of  the  neutral  wind.  In  general,  the  velocity  of  the  plasma  cloud  can  be 
written  in  the  form 


-♦  *-4 

E x B 
a 

+ TJ 

B 


E x B 
a 

B2 


(4. 16) 


where  the  coupling  parameter,  r\,  is  a function  of  x and  £,  the  FLIP  conductivity 

ratio  and  length  to  width  ratio,  respectively,  of  the  plasma  cloud.  The  functional 

dependence  of  r?  on  x and  0 is  most  easily  obtained  by  solving  Eq.  (4. 15)  where 

the  electric  field  has  the  constant  value  = 2 + V x § far  from  the  cloud. 

a a n 

The  functional  form  of  tj(x,  j3)  depends  on  the  density  profile  in  the  plasma 
cloud.  A simple  solution  of  Eq.  (4. 15)  exists  for  plasma  clouds  having  a uniform 
density  and  a sharp  edge  corresponding  to  a = 0 for  the  shape  defined  in  Section  2. 
When  the  plasma  cloud  boundary  has  an  elliptical  shape  with  major  axis  a 0 ^ in 
the  direction  of  the  neutral  wind  and  minor  axis  a 0 , the  electric  field,  E , 

inside  the  plasma  cloud  has  a constant  value. 


Ec  = 


E' 

a 


1 + 


0x 

TP? 


(4. 17) 


The  cloud  moves  with  the  constant  velocity  V = V + E * x B/B^.  By  making  use  of 

c n c 
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the  expression  (4. 17)  and  subsituting  into  Eq.  (4. 16),  we  derive  the 
expression  for  the  coupling  constant 


i?  (x,  0) 


1 + j3  + /3x  ' 


(4. 18) 


We  invoke  the  principle  that  the  tendency  to  striate  depends  on  the 
difference  in  velocity  between  a plasma  cloud  with  a finite  value  of  ft  and  a sheet- 
like cloud  which  corresponds  to  a value  of  f}  = • . In  particular,  we  claim  that  the 
time  for  an  elliptically-shaped  cloud  to  bifurcate,  striate,  or  develop  smaller 
scale  striations  is  proportional  to  the  time  it  takes  a sheet-like  configuration  of 
ionization  to  move  a characteristic  distance  in  the  direction  of  elongation  of  the 
cloud  with  respect  to  the  cloud.  We  will  call  this  time  a bifurcation  time,  t^,  and 
indicate  that  it  is  proportional  to 


= C 


i vc(»)  - vc(p) 


(1  + P + P x) 


]■ 


(4. 19) 


Thus,  the  time  to  develop  or  bifurcate  is  proportional  to  the  scale  size  of  the  stri 
ation  and  inversely  proportional  to  the  neutral  wind  velocity  in  the  frame  moving 
with  the  velocity  of  the  ambient  ionization.  Table  4. 2 presents  numerical  values 
of  the  term  in  brackets  for  representative  values  of  x and  0.  We  see  that  large 
values  of  x and  ff  slow  down  the  rate  of  development  of  plasma  clouds  or  large 
striations. 


Table  4. 2.  Coefficient  of  bifurcation  time. 
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The  zeroth-order  one-layer  model  equations  indicate  that  the  rod-like  tip 
of  a striation  should  behave  exactly  as  an  original  plasma  cloud  behaves,  proceed- 
ing through  the  same  stages  of  development.  Thus  one  would  expect  that  the  back- 
side of  the  tip  of  the  striation  would  steepen  and  then  it  should  bifurcate  and/or 
generate  a number  of  smaller  scale  sheets,  each  of  which  would  undergo  its  own 
similar  development.  The  time  scale  for  these  subsequent  developments,  being 
proportional  to  the  scale  size,  a,  of  the  cloud,  would  rapidly  diminish  and  in  a 
finite  time  all  striations  would  become  infinitely  thin.  The  scaling  indicated  by 
Eq.  (4. 19)  has  been  supported  by  observations  of  the  time  it  takes  for  large  barium 
plasma  clouds  to  striate.  The  striation  onset  time  for  cylindrical  clouds  tended  to 
vary  with  scale  size,  a,  neutral  slip  velocity  V , and  conductivity  ratio,  x, 
according  to  Eq.  (4. 10)  with  £ = 1 and  ~ 1.  5.  For  instance,  in  the  case  of  the 
Spruce  ion  cloud  with  a characteristic  radius  3 km,  Vg  = 50  m/s,  and  an  estimated 
conductivity  ratio  x ~ 3-4,  definite  structure  was  seen  in  the  ion  cloud  at  around 
11  minutes  after  release.  In  further  support  of  the  reduction  in  development  time 
with  scale  size,  by  22  minutes  following  the  release  of  the  cloud  (two  striation  on- 
set times),  the  striations  were  well -developed  with  many  tens  of  striations  having 
scale  sizes  in  the  hundreds  of  meters  range.  It  is  obvious  that  this  subsequent 
development  of  the  cloud  after  the  initial  formation  of  the  smaller  sizes  proceeded 
much  more  rapidly  than  the  development  of  the  original  cloud. 

However,  this  simple  description  ceased  to  apply  at  later  times  for  the 
smaller-scale  striations.  Many  striations  with  a scale  size  of  order  a ~ 200  m 
were  seen  to  maintain  their  identities  for  several  hundreds  of  seconds.  Further- 
more, the  STRESS  test  series  has  shown  that  many  fine-scale  striations  with  char- 
acteristic dimensions  of  order  hundreds  of  meters  can  exist  for  as  long  as  two 
hours  following  a release  and  thus  were  able  to  maintain  their  identities  for  many 
thousands  of  seconds  without  undergoing  further  development.  It  is  apparent  that 
at  the  small  scale  sizes,  the  zeroth-order  equations  are  inadequate  to  describe  the 
behavior  of  the  rod-like  striations.  They  apparently  become  stable  against  the 
formation  of  striations. 

There  are  a number  of  factors  that  have  not  been  included  so  far  in  the 
zeroth-order  model  that  could  be  important  for  affecting  the  development  of  small- 
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scale  striations.  Prominent  among  these  missing  factors  are  connections  to 
other  regions  of  the  ionosphere  along  the  same  magnetic  field  line,  the  proper 
treatment  of  diffusion  coefficients,  and  possibly  changes  in  the  effective  back- 
ground ionospheric  conductivity.  The  proper  treatment  of  the  mechanisms  that 
control  the  behavior  of  late-time  striations  may  also  lead  to  a proper  treatment 
of  the  factors  that  determines  the  initial  scale  size  in  larger  plasma  clouds. 


5.  bifurcation  model  and  time-dependent 

POWER-SPE  CTRAL  - DENSITIES 


5.  1 BIFURCATION  MODEL  APPLIED  TO  HANE 

In  Section  5. 2 below,  we  present  the  results  of  calculations  that  give 
the  phase  power -spectral -density  as  a function  of  time  based  on  parameters  asso- 
ciated with  HANE  and  the  concepts  described  in  Section  4.  In  this  section  we 
introduce  a methodology  for  describing  the  time -dependence  of  the  number  and 
scale  size  of  striations  that  might  be  expected  following  a HANE.  We  emphasize 
that  it  is  the  methodology  that  is  of  importance  rather  than  the  explicit  results 
produced.  However,  there  are  some  features  of  the  results  shown  in  this  and  the 
next  section  that  are  likely  to  occur  following  HANE. 

In  Section  4. 1 we  presented  a discussion  that  motivated  us  to  postulate 
that  around  500  s following  burst  there  exist  some  number,  Mq,  of  well-developed 
striations  with  scale  sizes  in  the  range  of  several  kilometers  to  several  tens  of 
kilometers.  In  this  section,  we  describe  a model  for  determining  the  time- 
dependence  of  the  number  and  scale  size  of  these  initial  striations  beyond  the 
first  500  s following  burst.  The  concept  is  simple.  For  any  striation  of  scale- 
size  a,  there  is  a characteristic  bifurcation  time  T^a)  • After  a time  lapse  of 
r ^(a),  we  envision  that  a large  striation  will  have  bifurcated  and  that  there  will 
then  be  two  striations  having  the  same  shape  and  peak  electron  density  as  their 
original  parent  striation.  In  order  to  conserve  the  number  of  electrons,  the  size 
of  these  daughter  striations  will  be  a /J~2  . 

We  now  define  a time -dependent  distribution  function  of  striaticn  scale 
sizes,  M(a,  t),  such  that  the  number  of  striations  of  scale  size  a within  the 
range  a to  a + da  is  M(a,  t)da.  If  we  label  the  initial  size  of  the  striation  a' 
and  the  present  size  a,  then  the  above  description  of  the  bifurcation  process  is 
consistent  with 


Substituting  a’  = JTa. , we  find  that  the  distribution  function,  M(a,  t)  would  obey  the 
equation 


M(a,  t)  = 2*M(,/Ta,  t - rb  (^a)  ) . (5.1) 

A specification  of  the  function  r^a)  is  necessary  in  order  to  find  a solution. 

Equation  (5. 1)  is  consistent  with  the  following  description.  A newly 
created  striationof  size  a'  exists  for  a time  T^a')  and  then  instantaneously  breaks 
into  two  striations  each  of  size  a’/ . This  description  of  the  bifurcation  process 
is  obviously  not  physical.  It  is  also  difficult  to  treat  analytically.  Instead,  we  will 
develop  another  equation  that  is  consistent  in  concept  with  the  above  bifurcation 
model  but  represents  striations  as  changing  continuously  in  size.  We  assume  that 
a striation  of  size  a after  a time  lapse  t was  earlier  part  of  a larger  striation 
of  size  a'.  Then  its  present  size  is  a function  of  the  earlier  size  and  the  time 
lapse  t,  i.  e. , a = a(a\  t).  A time  constant  T(a)  for  the  rate  of  change  of  size  of 
striation  of  size  a,  can  be  defined  by 
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1 _ 1 &a(a’.t) 

t (a)  a &t 


(5.2) 


where  we  take  t( a)  as  being  proportional  to  Tb(a)  . This  description  says  that  each 
striation  slowly  shrinks  continuously  creating  a fraction  of  a striation  of  the  same 
shrinking  size.  This  description  is  not  physical  either  but  it  roughly  satisfies 
Eq.  (5.  I),  it  is  analytically  tractable,  and,  for  a continuous  distribution  of  a large 
number  of  striations,  should  provide  a good  indication  of  the  si*> sequent  distribution 
of  scale  sizes  after  a number  of  bifurcations. 

If  all  striations  have  the  same  density  profile,  then  the  density  within  a 
striation  can  be  written  n(r)  = nQ  f (r/a).  In  order  that  the  bifurcation  process  con- 
serves the  number  of  electrons,  we  must  have 


K 


r 
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This  equation  determines  M(a,  t)  in  terms  of  the  initial  distribution  of  scale 
sizes  MQ(aT)  = M(a' , 0)  as 

M(a,  t)  = — Mo(a'(a,t)) 

a t 


/ a'(a,  t)\3  T(a)  „„  , , , , 

\ a ) t (av  (a,  t))  Mo(a(a,t))- 


(5.3) 


The  function  a'(a,t)  is  the  inverse  of  the  function  a(a',t)  and  represents  the 
initial  scale  size  of  the  striation  that  has  a size  a after  a time  lapse  t. 

Equation  (5.  3)  is  a general  result  for  the  time -dependent  striation  scale  - 
size  distribution  function  once  the  function  a’ (a,  t)  (or  equivalently,  a(ajt))  is 
determined.  In  order  to  provide  an  explicit  example  we  need  to  specify  an  appro- 
priate bifurcation  time  constant.  We  choose  to  illustrate  a simple  example  by 
considering  only  cylindrically-shaped  rod-like  striations  which  are  characterized 
by  a = 0 and  j3  = 1.  Motivated  by  the  discussion  following  Eq.  (4. 19),  we 
choose 


Tb(a)  = 


(a-asr 


(5.4) 


The  second  term  in  square  brackets  represents  an  ad  hoc  stabilizing  term  for 

small  scale  sizes  such  that  r^a)  - » as  a - a . We  shall  arbitrarily  choose 

the  minimum  scale  size  a = 0. 1 km  in  the  following  calculations.  For  moderate 

values  of  the  constant  C,  we  find  that  for  a ~ 3 km  and  V - 50  m/s,  the 

s 

bifurcation  time  given  by  Eq.  (5.  4)  is  ~ 10  min  which  is  about  the  time  that  it 
took  the  Spruce  ion  cloud  to  develop  its  initial  structure.  At  late  time  after 
many  striations  had  formed  in  the  Spruce  ion  cloud,  we  found  that  no  striation 
with  scale  size  <300  m was  seen  to  bifurcate  in  times  ~ 200  s.  We  choose 
C = 16  in  order  to  be  consistent  with  that  fact. 
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A number  of  other  functions  for  T^a),  ad  hoc  or  physically  motivated, 
can  be  examined.  For  instance,  the  stabilizing  term  could  have  the  form 


, v 10a 
Tbl  ^ ” V 


1 + 


C2as 

(5^) 


(5.5) 


For  comparison  purposes,  the  two  functions  Tb(a)  and  Tbi(a)  with  C = 16  and 
CL  = 5 are  shown  plotted  in  Fig.  17  with  V = 500  m/s.  Henceforth,  we  will  adopt 
the  time  constant  r(a)  associated  with  r^Ca)  solely  because  the  striations  do  not 
accumulate  at  the  smallest  scale  size  ag  as  rapidly  as  for  T^a). 

We  take  the  time  constant  defined  in  Eq.  (5.  2)  to  have  the  same  form  as 
rb(a),  i.e. , 


r (a) 


Toa 


Ca  ' 

1 4- 


(a  - a ) 
' s' 


(5.6) 


If  we  integrate  Eq.  (5.  2)  for  a bifurcation  time  rb(a)  and  neglect  the  time-dependence 
of  a on  the  left-hand  side,  we  obtain  by  using  Eqs.  (5. 4)  and  (5.  6), 


Tb(a) 
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ToV 
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- In 


a/y? 
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(5.  7) 


from  which 


20 

V In  2 
s 


s 


(5.8) 


with  V expressed  in  km/s.  As  a characteristic  value  of  the  neutral  slip  velocity 
s 

in  the  HANE  late-time-frame,  we  take  V = 0.  5 km/s.  Equation  (5.  8)  then  gives 

s 

t = 58  s as  the  value  on  which  our  numerical  calculations  below  are  based, 
o 

With  r (a)  of  the  form  given  by  Eq.  (5. 6),  the  time  it  takes  for  a stri- 
ation  of  size  a to  be  formed  from  a striation  of  size  a’  can  be  found  by  integrating 
Eq.  (5.2); 
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Bifurcation  time  as  a function  of  striation  scale  size,  a,  for  two  models  discussed 
in  the  text. 


t(a,af)  = - 


(5.9) 


da 


- a + Ca 

s 


2 


1 

a'-a 

s 


For  initial  scale  sizes  a‘  > 10  km,  C = 16,  and  a = 0. 1,  the  last  term  in 

s’ 

Eq.  (5.  9)  is  less  than  one  second  in  magnitude  and  we  will  neglect  it.  Henceforth 

we  will  measure  time  t in  seconds  after  burst  and  will  assume  that  at  t = t = 500  s 

o 

there  exist  striations  with  some  initial  distribution  in  scale  sizes.  In  terms  of  the 
normalized  time  lapse  F given  by 


t 


t(a,  ar)/r  Q 


t - 500 

“55 ’ 


(5.  10) 


we  may  use  Eq.  (5.9)  to  explicitly  approximate  the  functions  a' (a,  t)  and  a(a’  ,t)  : 


a'  = a + t - 


0. 16 

- 0.  1 


(5.  11) 


= 0.  05  + i-  £af  - t + J{ a'  - 0.1-t)2  + 0.64  ] . (5. 12) 


In  order  to  specify  the  striation  scale-size  distribution  as  a function  of 

time,  we  need  to  specify  the  distribution  of  sizes  at  the  initial  time,  t and  the 

time  dependence  of  the  number  of  new  striations  that  are  being  generated  at  each 

scale  size  other  than  by  the  bifurcation  process.  The  latter  represents  a time- 

dependent  source  of  striations  each  of  which  would  then  evolve  in  scale  size.  We 

will  illustrate  the  methodology  by  considering  only  the  initial- value  problem  and 

neglecting  an  additional  source  of  newly-created  striations.  However,  the  second 

initial  distribution  that  ve  discuss  below  is  motivated  by  the  distribution  of  sizes 

that  results  from  the  specific  choice  of  r (a).  The  first  example  of  an  initial  size 

distribution  is  a delta  function,  i.  e. , initially  there  are  striations  of  only  one  size 

present.  Based  on  the  discussion  in  Section  4. 1. 1,  we  shall  assume  that  at  500  s 

after  burst  there  are  Id  striations  of  scale  size  20  km.  Thus  M .(af)  = M 6(a'  -20). 

o olo 

By  applying  the  general  formula  given  in  Eq.  (5.  3),  we  find 
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Mx(a,  t) 


a'  aa’ 


6 (a'  - 20) 


400M 


5(a-a,(t)  ) 


(5.  13) 


where  a^(t)  is  the  solution  of  Eq.  (5. 12)  with  a'  =20.  The  second  form  for 
lVL(a,  t)  indicates  that  at  any  time  only  one  scale  size  is  present  and  the  number 

i 2 

of  striations  varies  as  1/a  . We  note  that  the  total  number  of  electrons  is  con- 
served (time-independent)  and,  according  to  the  equation  leading  to  the  derivation 

/CO 

x f (x)  dx  . 
o 

The  above  model  is  consistent  with  the  presumed  instantaneous  creation 

of  a number  of  striations  of  a given  size  at  500  s.  If,  instead,  striations  were 

created  continuously  throughout  the  initial  500  s time  interval,  the  striations  into 

2 2 

which  they  evolve  would  have  a 1/a  distribution.  An  initial  1/a  distribution 

2 

of  large  scale  sizes  tends  to  maintain  itself.  The  1/a  behavior  is  a property 
of  the  choice  of  the  form  of  r(a).  For  intermediate  values  of  a such  that 
a « a « a’,  we  find  from  Eq.  (5.  11)  that  a* (a,  t)  ~ t . By  applying  the 
general  formula  Eq.  (5.3),  we  obtain 


M(a,  t) 


= / a’  V r(a) 

\ a / 77a7) 


M (a') 
ov 


t2M  (t) 
o 


(5.  14) 


i.  e.,  M(a,  t)  <*  1/a  independent  of  the  form  of  the  initial  distribution  MQ(a’). 

A different  choice  of  the  function  t (a)  gives  a different  distribution.  For  example, 

— 3 

if  r (a)  = t,  a constant  independent  of  a,  then  M(a,  t)  « l/a  . 

2 

For  the  second  example  we  shall  take  an  initial  1/a  distribution. 

At  500  s after  release  we  assume  that  there  are  a total  of  M striations  distributed 

2 0 

in  size  between  10  and  40  kilometers  with  a 1/a  distribution.  Thus 


Mo2(a')  = 


H (a’  - 10)  H (40  - a') 


(5.  15) 
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where  H(x)  is  the  Heaviside  step-function.  The  numerical  coefficient  and  end  limits 
are  chosen  so  that  initially  there  are  the  same  number  of  striations  present  as  for 
M01(a'),  although  they  are  now  distributed  in  size,  and  they  contain  the  same 
number  of  electrons.  We  note  that  although  half  of  the  striations  present  have  scale 
sizes  between  10  and  16  km,  they  contain  only  20%  of  the  electrons  while  half  of 
the  electrons  are  contained  in  20%  of  the  striations  with  sizes  between  25  and  40  km. 

Thus  applying  the  general  formula  (Eq.  (5.3))  we  obtain  the  second  time- 
dependent  distribution  of  striation  scale  sizes  , 


Mg  (a,  t) 


40M 

o 


0. 16 
(a-0. 1): 


H(a-am(t))  H(aM(t)-a) 


(5. 16) 


o 

where  we  have  neglected  0. 16/(a’  - 0. 1)  compared  to  1 for  a'  > 10  km.  The  func- 
tions &m(t ) andaM(t  ) are  the  solutions  of  Eq.  (5. 12)  with  a'  = 10  and  40,  respec- 
tively. Note  that  Mg(a,  t)  has  a time -independent  (but  scale  -size  dependent)  ampli- 
tude but  has  a nonzero  value  only  in  the  finite  time -dependent  range  am(t)  < a <a^(t). 

Table  5. 1 lists  the  values  of  a.  (t),  a (t),  a,,(t),  and  the  number  of  striations, 

/oo  £ m m 

M.(a,  t)  da,  that  are  present  according  to  the  two  model  distributions  at 
- o ] 

500  s increments  after  burst . 


Table  5. 1 Bifurcation  model  parameters 
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a 
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a,  , 
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N1 

N2 

(s) 

(km) 

(km' 

(km) 

(M0) 

<N0) 

500 

20 

10 

40 

1 

1 

1000 

11.39 

1.49 

31.38 

3 

9 

1500 

2.82 

0.  122 

22.67 

50 

4,300 

2000 

0.  127 

0. 110 

14.  15 

24, 000 

13,  000 

00 

0.  100 

0.  100 

0.  100 

40,  000 

40,  000 

Figure  18  contains  a plot  of  the  amplitude  of  the  two  distributions  in 
scale  size  defined  by  Eqs.  (5. 13)  and  (5. 16)  respectively.  The  solid  curves  give 
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STR1ATI0N  SIZE  DISTRIBUTION  (M-/KM) 


STRIATION  SIZE  8 (KM) 

Figure  18.  The  two  distributions  in  striation  scale  size  resulting  from  the  two 
assumed  initial  distributions  corresponding  to  a delta  function. 

ancl  *°  a 1^a  distribution.  M o2(a).  as  discussed  in  the  text. 
* cXrt^ueS  * am’  aM'  and  al  as  defined  in  the  text  are  indicated  at 
t = 500  s,  1000  s,  1500  s,  and  2000  s after  burst,  respectively. 
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the  number  of  striations  of  a given  size.  For  M^(a,  t),  there  is  only  one  scale  size 
at  any  time.  The  value  of  that  scale  size  at  times  of  500,  1000,  1500  and  2000  s 
after  burst  are  indicated  by  small  vertical  bars  labeled  a^  through  a ^ , respec- 
tively. Likewise,  at  any  given  time  in  the  evolution  of  the  striations  according  to 

M0(a,  t)  there  are  always  striations  present  between  a minimum,  a , and  a 
z m 

1 maximum,  a..,  scale  size.  The  values  of  a and  a,,  at  500,  1000,  1500  and 

M’  m M 

2000  s after  burst  are  also  indicated  by  vertical  bars.  Figure  19  shows  the  total 
number  of  striations  present  according  to  the  two  models  as  a function  of  time 
t after  burst. 

« As  the  distribution  ML(a,  t)  evolves,  the  ratio  a,  ./a  at  first  increases 

[ and  then  decreases  as  can  be  seen  fromTable  5. 1.  At  2000  s after  burst  there  are 

still  large-scale  striations  with  sizes  greater  than  one  kilometer  present  as  well 
as  many  striations  close  to  the  minimum  size  a^.  After  a few  more  hundred 
seconds  aM  - ag  and  the  distribution  approaches  a delta  function  with  all  striations 
being  of  order  100  meters  in  size. 

5.2  TIME-DEPENDENT  PHASE  POWER-SPECTRAL-DENSITIES 

In  Section  3 of  this  report  it  was  emphasized  that  a knowledge  of  the 
phase  power-spectral-density  (PSD)  is  sufficient  for  determining  the  average  pro- 
perties of  scintillation  effects  to  be  expected  on  propagating  through  striated  media. 
In  Section  2,  the  effects  of  peak  electron  concentration,  scale  size  of  striations,  and 
shape  and  density  profile  of  striations  on  the  phase  PSD  produced  by  a single 
striation  was  calculated.  In  Section  4,  the  striation  characteristics  to  be  expected 
from  HANE  were  discussed  and  in  Section  5.  1 we  introduced  a simplified  bifurca- 
tion model  that  determined  the  number  of  striations  of  different  scale  sizes  as  a 
function  of  time.  In  this  section  we  compare  the  phase  PSDs  that  result  from  the 
two  initial  assumed  distributions  of  scale  sizes,  M^(a.  t)  and  MgCa.t).  We  will 
also  show  and  discuss  the  changes  in  the  PSD  as  the  size  distribution  of  striations 
changes  as  a function  of  time. 

The  results  shown  here  should  not  be  taken  as  a specific  example  of 
what  actually  occurs  in  HANE,  but  rather  as  a model  calculation  of  the  methodology 
that  we  have  been  describing.  However,  the  physical  characteristics  that  are 
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believed  to  operate  on  the  striations  during  HANE  will  undoubtedly  produce  some  of 
the  effects  that  will  be  shown  here,  As  we  learn  more  detailed  and  precise  informa- 
tion regarding  the  time  development  of  striations  in  HANE,  this  same  methodology 

can  be  applied  to  provide  more  accurate  estimates  of  the  time -dependent  PSD. 

o 

The  contribution  to  the  phase  PSD,  <p  (q),  due  to  a single  rod  (JJ  = 1) 
of  radius  a is  given  by  Eq.  (3. 1. 1)  with  a = 0.  The  PSD  resulting  from  a set  of 
rods,  located  randomly,  having  a distribution  of  scale  sizes  given  by  M^(a,  t)  is 

2 

(q,  t)  = C P.  (q,  t)  = 4f  ("T"2)  j a2  J^qa)  M . (a,  t)  da  (5. 17) 

where  P.(q,  t)  are  normalized  PSD.  By  expanding  the  Bessel  function  for  small 
q as  q -*  0,  we  define  the  constant  C to  have  the  value 

C a <p2  (0,500)  = 1.4x  105  »r  (rex  nQ)2  Mq/L 

2 

a (n  X)  M 9 

= 3. 5 x 10  j- rad  km  (5. 18) 

13  -2 

with  (nQX)  expressed  in  units  of  10  m and  L given  in  km.  With  this  choice, 
we  find  Pg  (0,  500)  = 1 and  P^  (0,  500)  = 4/7.  In  Figures  20  and  21  we  show  the 
time -development  of  P^(q,  t)  and  Pg  (q,  t),  respectively,  obtained  by  using  the 
two  different  scale-size  distribution  functions  M^(a,  t)  and  iV^a,  t)  discussed  in 
the  previous  section,  i.  e. , those  developing  from  an  initial  delta  function  and 
l/a  distribution,  respectively.  In  Figures  22through  25,  we  compare  P^(q,  t) 
with  Pg(q.  t)  at  times  corresponding  to  500,  1000,  1500  and  2000  s after  burst. 


Figure  20  shows  Pj(q,  t)  at  5 different  times  including  the  asymptotic 
curve  as  t -*  • . Because  P.(q,  t)  is  the  PSD  that  results  from  a collection 
of  400  M /a  (t)  striations  of  identical  scale  size  a , each  of  the  curves  shown 

O v v 

is  the  same  curve  just  displaced  vertically  and  horizontally  by  the  appropriate  { 

amount.  The  scales  on  the  bottom  show  the  appropriate  wavenumber,  q,  and 
the  equivalent  wavelength,  2 ir/q.  The  PSD  as  given  by  Eq.  (3. 1. 1)  has  zeros 
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Figure  20.  Time  development  of  the  Power -Spectral -Density  corresponding  to  an  initial  delta  function 
in  scale  size.  The  time  is  given  in  seconds  after  burst.  The  curve  labeled  500  s results 
from  the  assumed  initial  single  size  of  a - a (500)  ^ 20  km. 


at  the  zeros  of  the  Bessel  function.  The  first  five  zeros  are  indicated  only  for  the 
initial,  t = 500  s,  and  the  final  , t -»  ®,  curves.  For  the  other  three  intermediate 
times,  only  the  first  zero  is  indicated.  The  zeros  are  a manifestation  of  the  exis- 
tence of  a single  scale  size.  If  there  had  been  a slight  range  in  scale  sizes  about 
the  value  a , the  zeros  would  have  averaged  out  and  the  peaks  would  have  been 

reduced  in  amplitude.  The  higher  wavenumber  portions  of  the  PSD  have  been 

-3 

replaced  by  the  straight  line  representing  the  q asymptotic  curve  at  the  average 

values  of  the  function.  This  straight  line  is  located  a factor  of  2 below  the  peaks 

2 

because  the  average  value  of  sin  is  0.  5. 

The  values  of  P^(q,  t)  have  simple  expressions  in  the  limits  as  q -»  o 
and  q - ® . By  taking  the  appropriate  power  series  expansion  for  small  q and 
asymptotic  expansion  for  large  q of  the  Bessel  function,  we  find  that  P^q,  t)  can 
be  represented  fairly  accurately  by  a function  of  the  form 


Px(q,  t) 


' o 

rt  | 
i jl 

2 

[l- 

0,851  q2 

1 

(5.  19) 


2 4 

Because  P^(0,  t)  decreases  as  a^  (t),  its  value  at  t -*  ® is  a factor  of  4 x 10  less 

than  the  initial  value  at  t = 500  s.  The  phase  PSD  at  very  high  wavenumber  increases 

monotonically  as  1/a  (t)  amounting  to  a late -time  increase  of  a factor  of  200  over 

v 

the  initial  values  of  power  at  high  wavenumber.  From  the  expression  in  Eq.  (5.  19) 
one  finds  that  the  power  at  a given  wavenumber,  q.  is  a maximum  at  the  time  that 
a^t)  = 1.53/q. 


Figure  21  shows  a similar  plot  of  the  phase  PSD,  Pg(q.  t).  evaluated  at 
the  four  times  shown.  This  figure  shows  the  time  evolution  of  P2  resulting  from 
a group  of  striations  which  evolve  in  time  by  the  bifurcation  model  discussed  in  the 
previous  section,  having  an  initial  1/a  distribution  in  scale  sftses  ranging  from  a 
minimum  of  10  km  to  a maximum  of  40  km.  The  subsequent  figures  show  a 
direct  comparison  between  P^(q,  t)  and  P2(q,  t)  at  each  of  the  re  spective  times. 
The  first  observation  is  that  the  distribution  in  sizes  has  averaged  out  the  zeros  of 
the  Bessel  function  and  the  slight  ripple  seen  at  the  lower  wavenumbers  is  a 
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reflection  of  the  shape  of  the  size  distribution.  We  note  that  the  decrease  in  power 
at  low  wavenumbers  with  time  is  slower  for  Po(0,  t)  than  for  P.,(0,  t)  . This 

“ 1 2 

might  seem  surprising  at  first  since  P^O,  t)  is  initially  proportional  to  a (t) 

3 3^ 

while  the  initial  value  of  Po(0,  t)  is  proportional  to  a (t)  - a (t).  However, 

u ivl  m ^ ^ 

this  initial  small  decrease  is  easily  explained  by  realizing  first  that  a « a^  and 

q o nri  ivl 

initially  a^(t)  ~ (40  -1;  so  that  by  t = 2000  s,  (F  = 25.  9)  it  has  decreased  by 
only  a factor  of  less  than  25. 

At  times  after  1000  s,  Pgfait)  begins  to  take  on  an  interesting  character- 
istic. The  accumulation  of  a large  number  of  striations  around  a 100  meter  scale 
size  causes  an  enhancement  in  the  PSD  at  the  corresponding  characteristic  value 
of  q.  In  addition,  there  are  still  a large  number  of  electrons  contained  in  striations 
that  are  bigger  than  many  kilometers.  Thus,  it  appears  that  the  resultant  spectrum 
is  easily  explained  by  the  superposition  of  the  contributions  from  the  large  number 
of  electrons  contained  in  the  larger  striations  and  the  large  number  of  striations 
at  the  minimum  scale  size. 

Another  consequence  of  the  fact  that  a large  number  of  striations  are 
approaching  the  same  size,  approximately  100  m,  is  that  the  zeros  of  the  Bessel 
function  are  beginning  to  manifest  themselves  again.  Indeed,  in  the  limit  t -•  ®, 

Pgfq,  t)  approaches  P^(q,  t)  which  is  shown  as  the  fifth  curve  in  Figure  20.  The 
most  dramatic  change  that  occurs  in  the  phase  PSD  after  2000  s is  the  rapid 
decrease  in  power  in  the  low  wavenumbers. 

Figures  22  - 25  compare  the  two  phase  PSDs  from  the  two  distributions 
at  each  of  the  four  different  times.  Figure  22  is  a comparison  of  the  two  phase 
PSDs  resulting  from  the  assumed  initial  distributions.  Other  than  the  presence  of 
the  zeros  resulting  from  the  choice  of  a single  scale  size  only  for  M^,  it  is 
apparent  that  these  two  PSD  curves  are  comparable  to  each  other  to  within  a factor 
of  2.  At  1000  s shown  in  Figure  23,  the  differences  in  the  PSD  have  increased  in 
some  regions  of  q to  approximately  a factor  of  4 but  the  differences  betweem 
them  are  not  great  yet. 

Also  shown  in  Figure  23  is  a naive  application  of  the  effect  of  diffusion 
shown  by  a dashed  curve.  If  one  introduces  the  effect  of  diffusion  by  adding  a 
diffusion  term  to  the  continuity  equation  in  the  form  given  by  Eq.  (4.  7) , a 
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1000  s 


Hence,  both  the  linearized  and  nonlinear  analysis  applied  to  the  diffusion  equation 

with  a constant  diffusion  coefficient  results  in  the  fact  that  the  amplitude  of  the 

2 

Fourier  component  is  time  dependent  and  decreases  in  time  as  exp  (-q  Dt). 

If  this  indeed  were  the  effect  on  the  amplitude  of  the  Fourier  components  of  the 

2 

density  distribution,  the  phase  PSD  would  be  modified  by  a factor  exp  (-2  q Dt). 
The  dashed  line  shown  in  Figures  23-25  is  nothing  more  than  an  application  of  this 
factor  to  the  undisturbed  curves. 

We  emphasize  that  the  purpose  of  including  these  dashed  lines  is  to 
show  solely  what  a naive  application  of  diffusion  to  the  PSD  would  produce.  We  do 
not  feel  that  this  is  the  correct  method  of  showing  the  effect  of  diffusion  on  phase 
PSDs.  The  proper  method  of  including  diffusion  is  complicated  because  it  needs 
to  be  applied  to  a time-evolving  situation  and  the  effect  on  the  profile  of  striations 
needs  to  be  evaluated.  Part  of  the  reason  why  this  task  is  complicated  is  because 
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I 

the  value  of  diffusion  at  a particular  altitude  is  affected  by  the  electrodynamic  j 

coupling  to  other  regions  of  ionization  on  the  same  magnetic  field  line.  The  effects  i 

shown  by  the  dashed  line  grossly  overestimate  the  effect  that  diffusion  has  on  the  j 

phase  PSD  at  high  wavenumbers.  Due  to  the  presence  of  steepening  processes,  the  j; 

effect  of  diffusion,  properly  applied,  would  be  restricted  to  values  of  wavenumbers  j 

higher  than  the  value  qn  ~ 10  km  * indicated  in  these  figures.  For  example,  the  ; 

U I 

treatment  of  diffusion  in  Section  4. 1.  2 indicated  lengths  of  the  order  3 m (see  Eq.  i 

1 1 

(4. 1.  4) ) which  will  affect  only  values  of  q > 500  km  . \ 

Figure  24  shows  that  there  are  differences  in  the  phase  PSD  produced  by 
the  two  different  scale -size  distribution  functions  at  1500  s after  burst  that  approach  | 

an  order  of  magnitude  difference  in  power  for  some  values  of  q.  The  bifurcation  pro-  i 

i 

cess  described  in  the  previous  subsection  is  an  incompressible  process.  As  such,  !• 

i 2 

the  Jr  n (r)  dr  is  conserved  as  the  striations  bifurcate  and  change  their  number 

® A 2 ^ 2 j, 

and  size.  This  fact  implies  that  fn  (k  ,k  )dk  dk  is  conserved  where  n (k  ,k  ) I 

* J e v x*  y'  x y e ' x’  y'  i 

i is  the  PSD  of  the  fluctuations  in  electron  density.  For  an  isotropic,  homogeneous 

distribution  of  cylindrically- symmetric  rods,  the  case  that  we  have  been  considering, 

2*2  a 2 2 2 4 * 

the  phase  PSD  is  given  by  <o  (q)  « n (q.o),  n is  a function  of  k = (k  +k  ) , 

0 0 x y j 

and  the  conserved  quantity  in  k-space  is  J*  kne^(k)dk.  Thus,  the  value  of  the 

J qP.(q,  ^ ^q  f°r  the  two  curves  in  Figure  24  is  also  conserved.  While  the 

value  of  this  integral  over  all  q is  the  same  for  both  P^  and  Pg,  the  contribution 

to  the  integral  from  the  middle  range  of  q values  is  greater  for  P^  than  for  | 

Pg  has  relatively  larger  contributions  to  its  total  value  from  low  and  high  values  j 

of  q.  Note  that  J qP(q,  t)dt  is  not,  in  general,  conserved  for  asymmetric  striations 


w ith  0 > 1. 

Figure  25  shows  that  the  Pgfa,  t)  is  approaching  that  of  P j(q,  t)  as  time 
approaches  infinity.  The  large  power  at  low  values  of  q for  P2  at  this  time  is 
just  a reflection  of  the  fact  that  almost  half  of  the  electrons  still  reside  in  striations 
of  greater  than  1 kilometer  radius. 

5. 3 SUMMARY 

The  time -dependent  PSD  associated  with  a given  set  of  striations  is  suffi- 
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of  the  shape  of  striations  and  the  density  profile  within  striations  in  determining 
the  phase  PSD  have  been  incorporated  and  the  results  have  been  studied  in  terms 
of  the  profile  and  shape  parameters  a and  £ . 


The  weak  link  in  the  procedure  for  determining  the  propagation  effects 
that  result  from  HANE  is  in  being  able  to  specify  accurately  the  characteristics  of 
striations.  Based  on  our  currently  available  understanding  of  the  processes  which 
operate  to  create  the  structured  ionization,  we  have  presented  a discussion  of 
the  range  of  expected  striation  scale  sizes  that  can  be  expected  to  exist  500  s 
after  a high  altitude  nuclear  burst.  Furthermore,  we  have  discussed  our  current 
understanding  of  the  time  development  of  the  striations. 

A perhaps  surprising  result  is  that  at  late  times,  after  500  or  1000  s 
following  burst,  the  region  of  enhanced  ionization  has  broken  up  into  individual 
large  size  striations  whose  field-line-integrated  Pedersen  (FLIP)  conductivity  is 
of  the  same  order  of  magnitude  as  the  FLIP  conductivity  created  by  the  deploy- 
ment of  48  kg  barium  releases.  This  similarity  strongly  suggests  that  the  con- 
tinued development  of  the  striations  in  HANE  beyond  500  s after  burst  will  be 
similar  to  the  development  of  striations  in  barium  plasma  clouds  with  appropriate 
adjustments  in  the  time  scale  due  to  the  different  initial  size  and  neutral  slip 
velocity.  As  the  striations  break  up  into  smaller  scale  sizes  the  time  for  them 
to  continue  to  evolve,  or  bifurcate,  continually  decreases.  On  a time  scale 
measured  in  several  500  s units,  we  reach  the  conclusion  that  most  of  the  ioni- 
zation will  exist  in  striations  of  less  than  a kilometer  in  size.  Observations  of 
barium  clouds  suggest  that  there  is  a minimum  size  striation  that  is  produced  that 
does  not  subsequently  striate.  This  size  appears  to  be  of  order  100  meters  but 
may  be  somewhat  less.  We  do  not  presently  have  adequate  understanding  for 
determining  this  size. 

Consistent  with  the  above  description,  we  have  introduced  a bifurcation 
model  that  contains  some  of  the  physical  features  operating  to  change  the  scale 
size  of  striations.  We  have  accomplished  this  by  introducing  a model  for  the 
bifurcation  time  as  a function  of  striation  scale  size.  The  specific  model  that 
we  use  is  motivated  by  our  physical  understanding  for  large  scale  sizes  but  its 
dependence  near  the  smallest  scale  size  is  purely  ad  hoc.  However,  it  reproduces 
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the  features  of  limiting  the  further  bifurcation  of  the  small  size  striations  and 
appears  to  give  times  that  are  consistent  with  the  observations  of  barium  clouds 
for  striations  of  around  300  m radius. 


An  effect  of  having  an  initial  distribution  in  scale  sizes  is  that  the  PSD 
does  not  maintain  a simple  shape  as  the  striations  evolve.  The  methodology  that 
we  have  described  allows  us  to  introduce  appropriate  modifications  as  our  under- 
standing increases.  Once  we  are  more  certain  of  the  initial  distribution  in  scale 
sizes,  the  values  of  the  parameters  a and  p as  a function  of  time,  and  the 
mechanism  that  limits  the  smaller  scale  sizes,  this  knowledge  can  be  incorporated 
into  the  procedures  that  we  have  discussed. 

Based  on  an  initial  exploration  of  the  appropriate  value  ranges  of  the 
parameters,  we  summarize  the  propagation  effects  in  HANE.  The  ellipticity  of  the 
striations  reduce  the  scintillation  effects  over  most  of  a striated  region  and  enhance 
the  effects  in  the  central  region.  Striations  with  high  central  electron  density  com- 
pared to  background  have  predominantly  parabolic  profiles  which  reduces  the  scin- 
tillation effects  to  that  of  abrupt  profiles  with  lower  electron  densities  closer  to 
the  background  in  most  cases.  High  frequencies,  1 s3  cm,  may  have  negligible 
scintillation  many  minutes  after  burst.  Lower  frequencies  may  have  Rayleigh 
scintillation  for  hours  after  burst. 

In  conclusion,  the  areas  in  which  additional  work  is  needed  are:  an 
understanding  of  the  physics  that  determines  the  initial  distribution  of  scale  sizes 
in  a cylindrical  geometry;  a more  accurate  understanding  of  the  time -development 
of  the  striations  in  the  nonlinear  phase;  the  mechanisms  that  lead  to  the  cessation 
of  the  bifurcation  process,  and  the  influence  that  the  steepening  and  diffusion  pro- 
cesses have  on  determining  the  density  profile  of  the  striations.  This  latter 
investigation  is  important  for  determining  the  wavenumbers  at  which  the  PSD 
decreases  drastically. 
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ATTN:  Code  810,  J.  Barna 
ATTN:  Code  480 
ATTN:  Code  101B 

ATTN:  Code  R1033,  M.  Raffensperger 

Defense  Documentation  Center 
12  cy  ATTN:  DO 

Defense  Intelligence  Aqency 
ATTN:  0C-7D,  W.  Wittig 
ATTN:  DB-4C,  E.  O’Farrell 
ATTN:  HQ-TR,  J.  Stewart 
ATTN:  DB,  A.  Wise 
ATTN:  DT-5 

ATTN:  DT-1BZ,  R.  Morton 
ATTN:  DT-1B 


Defense  Nuclear  Agency 


ATTN: 

DOST 

ATTN: 

STVL 

3 cy  ATTN: 

RAAE 

4 cy  ATTN: 

TITL 

Field  Command, 

Defense  Nuclear  Aqency 

ATTN: 

FCPR 

Interservice 

Nuclear  Weapons  School 

ATTN: 

TTV 

Joint  Chiefs 

of  Staff 

ATTN: 

J-3,  WWMCCS  Evaluation 

Joint  Strat. 

Tgt.  Planninq  Staff 

ATTN: 

JPS1,  G.  Goetz 

ATTN: 

JLTW-2 

ERADCOM  Technical  Support  Activity 
Department  of  the  Army 

ATTN:  DELET-ER,  H.  Bomke 

Harry  Diamond  Laboratories 
Department  of  the  Army 
ATTN:  DELHD-N-NP 
ATTN:  DELHD-N-RB , R.  Williams 
ATTN:  DELHD-N-TI , M.  Weiner 

U.S.  Army  Ballistic  Research  Labs. 

ATTN:  Tech.  Library, 

U.  S.  Army  Comm-Elec.  Engrq.  Instal . Agency 
ATTN:  CCC-EMEO; 

U.S.  Army  Foreign  Science  8 Tech.  Center 
ATTN:  DRXST-SD 

U.S.  Army  Materiel  Dev.  8 Readiness  Command 
ATTN:  DRCLDC,  J.  Bender 

U.S.  Army  Nuclear  & Chemical  Agency 
ATTN:  Library 

U.S.  Army  Satellite  Comm.  Agency 
ATTN:  Document  Control 

U.S.  Army  TRADOC  Systems  Analysis  Activity 
ATTN:  ATAA-PL 

ATTN:  ATAA-TCC , F.  Payan,  Jr. 

ATTN:  ATAA-TDC,  J.  Hesse 

DEPARTMENT  OF  THE  NAVY 

Naval  Electronic  Systems  Command 

ATTN:  NAVELEX  3101,  T.  Huohes 
ATTN:  PME  117 
ATTN:  Code  5011 
ATTN:  PME  117-T 
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DEPARTMENT  OF  the  navy  (fnnfin.IOH) 


Nava!  Intelligence  Support  Center 
ATTN:  N I SC-50 

Naval  Ocean  Systems  Center 

ATTN:  Code  8151,  C.  Baqqett 
3 cy  ATTN:  Code  5324,  W.  Moler 

Naval  Research  Laboratory 
ATTN:  Code  7580 
ATTN:  Code  7555 
ATTN:  Code  6700,  T.  Coffey 
ATTN:  Code  6707,  J.  Davis 
?tT™:  Code  7500,  Hg.  Corrni.  Dir.,  B.  Wald 
ATTN.  Code  6701,  J.  Brown 

Naval  Space  Surveillance  System 
ATTN:  J.  Burton 

Naval  Surface  Weapons  Center 

ATTN:  Code  F14,  R.  Butler 

Office  of  Naval  Research 
ATTN:  Code  420 
ATTN:  Code  421 

0ffiCeA?™-he0pC-60e4f  °f  NaVa1  °Peratl0nS 
ATTN:  Op-941 

Strategic  Systems  Project  Office 

Department  of  the  Navy 

ATTN:  NSSP-2722,  F.  Wimberly 
ATTN:  NSP-2141 

Naval  Surface  Weapons  Center 

White  Oak  Laboratory 
ATTN:  Code  F31 

DEPARTMENT  OF  THE  AIR  FQRCF 

Aerospace  Defense  Command/DC 

Department  of  the  Air  Force 
ATTN:  DC,  Mr.  Long 

Aerospace  Defense  Command/XPD 

Department  of  the  Air  Force 
ATTN:  XPDQ 
ATTN:  XP 

Air  Force  Avionics  Laboratory 
ATTN:  AAD,  W.  Hunt 
ATTN:  AAD,  A.  Johnson 

Air  Force  Geophysics  Laboratory 
ATTN:  LKB,  K.  Champion 
ATTN:  0PR1,  J.  Ulwick 
ATTN:  PHP,  J.  Aarons 
ATTN:  PHI,  J.  Buchau 
ATTN:  PHD,  J.  Mullen 
ATTN:  OPR,  A.  Stair 

Air  Force  Weapons  Laboratory 
ATTN:  DYC,  J.  Barry 
ATTN:  SUL 
ATTN:  DYC,  J.  Kama 

Deputy  Chief  of  Staff 

Research,  Development,  S Acq. 

Department  of  the  Air  Force 
ATTN:  AFRDQ 


3i£M01ENT_OF  THE  AIR  FORCE  (Continued! 

Electronic  Systems  Division,  AFSC 
ATTN:  DCKC,  J.  Clark 
ATTN:  XRW,  J.  Deas 
ATTN:  YSEA 

Foreign  Technology  Division,  AFSC 
ATTN:  TQTD,  B.  Ballard 
ATTN:  NIIS,  Library 

Rome  Air  Development  Center,  AFSC 
ATTN:  V.  Coyne,  OCSE 
ATTN:  Documents  Library/TSLD 

Rome  Air  Development  Center,  AFSC 
ATTN:  EEP  - 

Space  & Missile  Systems  Organization/MN 
Air  Force  Systems  Conroand 

ATTN:  MNNL,  S.  Kennedy 

Space  A Missile  Systems  Organization/SK 
Air  Force  Systems  Coronand 
ATTN:  SKA,  M.  Clavin 

Strategic  Air  Command/XPFS 
Department  of  the  Air  Force 
ATTN:  NRT 

ATTN:  DCX,  Chief  Scientist 

DEPARTMENT  OF  FNFPGV 

Department  of  Energy 
Albuquerque  Operations  Office 

ATTN:  Doc.  Con.  for  D.  Sherwood 

Department  of  Energy 
Library  Room  G-042 

ATTN:  Doc.  Con.  for  A.  Labowitz 

Office  of  Military  Application 
Department  of  Energy 

ATTN:  Doc.  Con.  for  D.  Gale 

OTHER  GOVERNMENT  AGFNf IFS 

Central  Intelligence  Agency 

ATTN:  RD/SI,  Rm.  5G48,  Hq.  Bldg 
for  OSI/PSTD 

Department  of  Commerce 
National  Bureau  of  Standards 

ATTN:  Sec.  Officer  for  R.  Moore 

Department  of  Transportation 
Office  of  the  Secretary 
ATTN:  R.  Lewis 
ATTN:  R.  Doherty 

Institute  for  Telecommunications  Sciences 
National  Teleconmunications  & Info.  Admin. 

ATTN:  L.  Berry 
ATTN:  A.  Jean 
ATTN:  W.  Utlaut 
ATTN:  D.  Crombie 

National  Oceanic  i Atmospheric  Admin. 
Environmental  Research  Laboratories 
Department  of  Coronerce 

ATTN:  Aeronoroy  Lab.,  G.  Reid 
ATTN:  R.  Grubb 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Aerospace  Corp. 

General  Electric  Company- TEMPO 

ATTN:  F.  Morse 

Center  for  Advanced  Studies 

ATTN:  N.  Stockwell 

ATTN:  W.  Knapp 

ATTN:  I.  Garfunkel 

ATTN:  M.  Stanton 

ATTN:  D.  01 sen 

ATTN:  B.  Gambill 

ATTN:  S.  Bower 

ATTN:  T.  Stevens 

ATTN:  J.  Carter 

ATTN:  D.  Chandler 

ATTN:  T.  Salmi 

ATTN:  DAS I AC 

ATTN:  R.  Slaughter 

ATTN:  V.  Josephson 

General  Electric  Tech.  Services  Co., 
HMES 

Inc. 

Analytical  Systems  Engineering  Corp. 

ATTN:  G.  Millman 

ATTN:  Radio  Sciences 

General  Research  Corp. 

Berkeley  Research  Associates,  Inc. 

Santa  Barbara  Division 

ATTN:  J.  Workman 

ATTN:  J.  Ise,  Jr. 

ATTN:  J.  Garbarino 

Boeing  Co. 

ATTN:  J.  Kenney 

Geophysical  Institute 

ATTN:  G.  Keister 

University  of  Alaska 

ATTN:  D.  Murray 

ATTN:  N.  Brown 

ATTN:  G.  Hall 

ATTN:  T.  Davis 

ATTN:  S.  Tashird 

ATTN:  Technical  Library 

Brown  Engineering  Company,  Inc. 

GTE  Syl vania,  Inc. 

ATTN:  R.  Deliberis 

Electronics  Systems  Grp. -Eastern  Div 
ATTN:  M.  Cross 

University  of  California  at  San  Diego 

ATTN:  H.  Booker 

HSS,  Inc. 

ATTN:  D.  Hansen 

Charles  Stark  Draper  Lab.,  Inc. 

ATTN:  D.  Cox 

University  of  Illinois 

ATTN:  J.  Gilmore 

ATTN:  K.  Yeh 

Computer  Sciences  Corp. 

Institute  for  Defense  Analyses 

ATTN:  H.  Blank 

ATTN:  E.  Bauer 

ATTN:  J.  Aein 

Comsat  Labs. 

ATTN:  J.  Bengston 

ATTN:  G.  Hyde 

ATTN:  H.  Wolfhard 

ATTN:  R.  Taur 

International  Tel.  & Telegraph  Corp. 

Cornell  University 

ATTN:  Technical  Library 

Department  of  Electrical  Engineering 

ATTN:  D.  Farley,  Jr. 

JAY COR 

ATTN:  S.  Goldman 

Electrospace  Systems,  Inc. 

ATTN:  H.  Logston 

Johns  Hopkins  University 

Applied  Physics  Lab. 

ESI,  Inc. 

ATTN:  T.  Potemra 

ATTN:  C.  Prettie 

ATTN:  Document  Librarian 

ATTN:  J.  Roberts 

ATTN:  J.  Marshall 

Kaman  Sciences  Corp. 

ATTN:  T.  Meagher 

Ford  Aerospace  & Comnunications  Corp. 

ATTN:  J.  Mattingley 

Lawrence  Livermore  Laboratory 
University  of  California 

General  Electric  Co. 

ATTN:  Doc.  Con.  for  Technical 

Informati 

Space  Division 

Dept.  Library 

ATTN:  M.  Bortner 

Linkabit  Corp. 

General  Electric  Co. 

ATTN:  I.  Jacobs 

ATTN:  F.  Reibert 

Lockheed  Missiles  & Space  Co.,  Inc. 
ATTN:  Dept.  60-12 

ATTN:  D.  Churchill 
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DEPARTMENT 

OF  DEFENSE  CONTRACTORS  (Continued 

Lockheed  Missiles  and  Space  Co.,  Inc. 

Rand  Corp. 

ATTN 

R.  Johnson 

ATTN 

E.  Bedrozian 

ATTN 

M.  Walt 

ATTN 

C.  Crain 

ATTN 

W.  Imhof 

Riverside  Research  Institute 

Los  Alamos 

Scientific  Laboratory 

ATTN 

V.  Trapani 

ATTN 

Doc.  Con.  for  D.  Westervelt 

ATTN 

Doc.  Con.  for  R.  Taschek 

Sandia  Laboratories 

ATTN 

P.  Keaton 

Livermore  Laboratory 

ATTN 

Doc.  Con.  for 

T.  Cook 

M.I.T.  Lincoln  Lab. 

ATTN 

Doc.  Con.  for 

B.  Murphey 

ATTN 

D.  Towle 

ATTN 

L.  Loughlin 

Sandia  Laboratories 

ATTN 

D.  Clark 

ATTN 

Doc.  Con.  for 

3141 

ATTN 

P.  Waldron 

ATTN 

Doc.  Con.  for 

Space  Project  Div 

ATTN 

Doc.  Con.  for 

D.  Thornbrouah 

McDonnell  Douglas  Corp. 

ATTN 

Doc.  Con.  for 

W.  Brown 

ATTN 

N.  Harris 

ATTN 

Doc.  Con.  for 

D.  Dahlgren 

ATTN 

G.  Mroz 

ATTN 

J.  Moule 

Science  Applications,  Inc. 

ATTN 

W.  Olson 

ATTN 

L.  Linson 

ATTN 

D.  Hamlin 

Mission  Research  Corp. 

ATTN 

D.  Sachs 

ATTN 

R.  Hendrick 

ATTN 

E.  Straker 

ATTN 

S.  Gutsche 

ATTN 

J.  McDougall 

ATTN 

F.  Fajen 

ATTN 

C.  Smith 

ATTN 

D.  Sowle 

ATTN 

D.  Baxter 

ATTN 

R.  Bogusch 

ATTN 

P.  Fischer 

Science  Applications,  Inc. 

ATTN 

M.  Scheibe 

Huntsvi 1 le 

Division 

ATTN 

D.  Divis 

Mitre  Corp. 

ATTN 

C.  Callahan 

Science  Applications,  Inc. 

ATTN:  G.  Harding 

ATTN 

SZ 

Mitre  Corp 

SRI  International 

ATTN 

W.  Foster 

ATTN 

G.  Carpenter 

ATTN 

M.  Horrocks 

ATTN 

D.  Johnson 

ATTN 

W.  Hall 

ATTN 

C.  Rino 

ATTN 

G.  Smith 

Pacific-Sierra  Research  Corp. 

ATTN 

G.  Price 

ATTN 

E.  Field,  Jr. 

ATTN 

R.  Leadabrand 

ATTN 

J.  Peterson 

Pennsylvania  State  University 

ATTN 

W.  Jaye 

Ionosphere 

Research  Lab. 

ATTN 

W.  Chesnut 

ATTN 

Ionospheric  Research  Lab. 

ATTN 

M.  Baron 

ATTN 

A.  Burns 

Photometries,  Inc. 

ATTN 

D.  Neil  son 

ATTN 

I.  Kofsky 

Physical  Dynamics,  Inc. 
ATTN:  A.  Thompson 


Tri-Com,  Inc. 

ATTN:  D.  Murray 


Physical  Dynamics,  Inc. 
ATTN:  E.  Tremouw 


TRW  Defense  & Space  Sys.  Group 


RAD  Associates 


ATTN 

ATTN 

ATTN 


S.  Altschuler 
R.  Plebuch 
D.  Dee 


ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 

ATTN 


R. 

E. 

R. 

F. 
H. 
W. 
W. 

c. 


Lelevier 
Gabbard 
Turco 
Gi  lmore 
Ory 

Karzas 
Wright,  Jr. 
MacDonald 


Vi  si  dyne,  Inc. 

ATTN:  J.  Carpenter 
ATTN:  C.  Humphrey 
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